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Abstract 

We generalize Ringstrom's global future causal stability results 
|14j for certain expanding cosmological solutions of the Einstein-scalar 
field equations to solutions of the Einstein-Maxwell-scalar field sys- 
tem. In particular, after noting that the power law inflationary space- 
times {M^~^^ , g , (j)) considered by Ringstrom in [13] are solutions of 
the Einstein-Maxwell-scalar field system (with exponential potential) 
as well as of the Einstein-scalar field system (with the same exponen- 
tial potential), we consider (nonlinear) perturbations of initial data 
sets of these spacetimes which include electromagnetic perturbations 
as well as gravitational and scalar perturbations. We show that if (as 
in [Hj) we focus on pairs of relatively scaled open sets Ujig C Uijig 
on an initial slice of (M""*"^,^), and if we choose a set of perturbed 
data which on U^Rg is sufficiently close to that of {M"'~^^ , g, ^, A = 
0), then in the maximal globally hyperbolic spacetime development 
{M^^^,g, (j), A) of this data via the Einstein-Maxwell-scalar field equa- 
tions, all causal geodesies emanating from [//jp are future complete 
(just as in {M^~^^,g)). We also verify the controlled future asymptotic 
behavior of the fields in the spacetime developments of the perturbed 
data sets. 
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1 Introduction 



1.1 Background 

The well-posedness of the Einstein vacuum equations, the Einstein-Maxwell 
equations, and other related field equation systems has been established for 
many years [U [2] • Specifically, it is known that for any given set of initial 
data satisfying the constraint equations, there exists a unique solution to the 
Einstein (or Einstein-Maxwell, etc.) equations for some amount of time to 
the future. Well-posedness also establishes that small perturbations to an 
initial data set only lead to small changes to the corresponding solution in 
finite time, and that if those changes to the data are confined to a subset 
of the initial hypersurface then the changes in the solution occur strictly in 
the domain of dependence of that subset. However, well-posedness gives no 
information about global- in-time behavior of the development. 

One way to formulate this issue for Einstein's and other relativistic equa- 
tions is in terms of global future causal stability ( "GFC-stability" ) which 
addresses the following question: Given a particular set of initial data for 
which the maximal globally hyperbolic development ( "MGH-development" ) 
is future causally geodesically complete, if one makes small perturbations to 
that initial data set, is the resulting MGH-development also future causally 
geodesically complete? Note that while the property of well-posedness gen- 
erally characterizes a PDE system together with all (or none) of its solutions, 
GFC-stability pertains to a particular solution or family of solutions of the 
system. Note also that GFC-stability concerns the full nonlinear PDE sys- 
tem (Einstein, Einstein-scalar, etc.), not a linearization of the system, and 
in referring to "perturbations to the initial data set" , we mean new data sets 
suitably near the original given set of initial data (which, for convenience 
here, we will refer to as the "background" data set). 

Global stability (GFC-stability or other related varieties) has been studied 
extensively for Einstein's theory, but has been established for only a very 
small number of solutions: The epic work of Christodoulou and Klainerman 
[6j (see also the later generalizations and simplifications in [31 [11]) proves the 
global stability of Minkowski spacetime for the vacuum and Einstein-Maxwell 
equations, while that of Friedrich [7| shows that the DeSitter spacetime is 
globally stable for Einstein's equations with a cosmo logical constant. More 
recently, Andersson and Moncrief [Ij have proven that the Milne spacetime^ 

"'^ These expanding spacetimes are constructed by spatially compactifying the mass hy- 
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are globally stable solutions of the vacuum equations. 

One feature of the DeSitter and the Milne solutions which makes it a bit 
easier to establish global stability for them is the fact that they are expanding 
solutions. In a rough sense, this property acts to inhibit the concentration of 
curvature, so that perturbations do not tend to lead to singularities forming. 
Hence in searching for solutions expected to be globally stable, one is led to 
consider expanding solutions. 

Motivated both by this consideration and by the recent astrophysical 
evidence |8] that our universe is likely expanding at an accelerated rate, 
Ringstrom has recently shown that certain solutions of the Einstein-scalar 
field equations with accelerating expansion are GFC-stable. He does this for 
both exponentially expanding background spacetimes satisfying the Einstein- 
scalar field equations with fairly general scalar field potential functions V{(j)) 
[13\ . and for power law expanding background spacetimes satisfying the 
Einstein-scalar field equations with a certain set of exponentially-decaying 
scalar potential functions [E]. 

In this work, we show that the power law expanding solutions considered 
by Ringstrom in [T?] are globally stable with respect to the Einstein- MaxweW- 
scalar field equations. By choosing the electromagnetic fields to vanish, we 
may consider Ringstrom's solutions from [l^ to be solutions of the Einstein- 
Maxwell-scalar system. To prove stability of these background solutions in 
the larger PDE system of course requires us to allow the perturbation solu- 
tions to include non-vanishing electromagnetic fields. This paper shows that 
this can be done, and that GFC-stability holds. We note that in a recent 
paper [15], stability has been proven for an electromagnetic generalization of 
Ringstrom's results for exponentially expanding solutions. 

A key feature of background solutions with sufficiently accelerated expan- 
sion is that the analysis can be strongly localized. This is because the entire 
future of a small subset in the initial hypersurface is determined completely 
by the initial data on a small neighborhood of that subset. Effectively then, 
the topology of the Cauchy slices of the background solutions being tested 
for stability and the topology of the perturbed solutions is irrelevant. 

The general structure of our proof is very similar to that of [H]: i) lo- 
calizing the analysis to the development of data sets on open sets in the 
initial Cauchy slice, with the formal extension of such local data sets to spa- 
tial tori; (ii) establishing the well-posedness of the Cauchy problem for the 

perboloids in Minkowski spacetime 
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field perturbations relative to the background solutions (with appropriate 
handling of the gauge choice); iii) defining energy-type functionals for the 
perturbation fields and their derivatives, and (with the help of bootstrap 
assumptions) proving monotonicity estimates for them; iii) using the energy 
estimates together with bootstrap arguments to prove long-time existence, 
regularity, and global estimates for the MGH-developments of the perturbed 
initial data; (iv) using the global estimates to analyze the dynamics of causal 
paths in the perturbed spacetimes, and thereby verifying future causal com- 
pleteness. The theme of this paper is showing that all of these steps work 
for the Einstein-Maxwell-scalar field equations (with exponentially-decaying 
scalar field potentials). 

1.2 Field Equations and Background Solutions 

Before stating our main results and proceeding to prove them, we wish to 
set up the field equations for the parametrized set of Einstein-Maxwell-scalar 
field theories which we work with here, and we wish to also state what the 
background solutions are, explicitly. 

The field variables for the Einstein-Maxwell-scalar field theories include 
the spacetime metric the electromagnetic vector potential A, and the scalar 
field 0. Letting R^y and R denote the Ricci tensor and the scalar curvature 
for g, letting F denote the electromagnetic tensor for A, and choosing the 
scalar field potential to take the form V{(j)) = Vqc"'^'^ (for constants Vq and 
A), we can write the field equations for this theory (for (n + l)-dimensional 
spacetimes) in the following (index) form: 

R^iu '^Rg^u 'l^fii/i (1) 

V'^V^0-V'(0) = O, (2) 
V^F,, = 0. (3) 

Here the stress-energy tensor for this system is given by 

T,. = d,(Pd,(P - g^,(}-g'"'dpci>d^ct> + V(0)) + {F,,F, - ^g^^F^^F^, (4) 

and V'{(f)) = —XV{(j)). Note that ([T]) can be rewritten as 

R^.u = 5^05,0 + ^--r (0)(?^, + F^^F, - ^ g^^F.^FP''. (5) 
n—1 / m — 1 
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Note also that three parameters characterize these systems of field equations: 
the spatial dimension n > 3, the scalar potential scale Vq > and the scalar 
potential decay coefficient A > 0; hence for convenience, we shall denote a 
particular choice of these theories by "Einstein-Maxwell-scalar{„ vb a}" • Fi- 
nally, note that the form of the field equations (II])- (H]) is consistent with the 
assumption here that the scalar fields are not charged, and so the interaction 
between the electromagnetic and the scalar fields is indirect (through the 
gravitational fields). 

We now wish to specify the background fields {g,(f),A), which i) are so- 
lutions of the system (II])- ii) have accelerating expansion, and iii) (as 
we shall show) are GFC-stable. The fields are defined on the manifold 
j\^n+i _ rpn ^ j^^^ Qj-^ which we choose the time coordinate t > and 
the global periodic spatial coordinates x*. If we now choose the constant 
parameters to > 0,p > 1, cq, and n, we write the following: 



-dt^ + e^^{t/toyP6ijdx'dx^, 

X 2, , Co 
0=-lnt--, 

A,, = 0. 



(6) 
(7) 
(8) 



These fields do not generally satisfy the field equations ([I])-®. However, 
if we require the field equation parameters {n, Vq, A} and the field parameters 
{to,p, Co, k} to satisfy the constraining relations 



A 



Co 



[(n- l)p]i/2 



In 



{n — l){np — l)p 
2Vo 



(9) 
(10) 



then indeed the fields [g, (j), A) do constitute a solution. Note that for a 
fixed spatial dimension n, ([9]) expresses a one-to-one correspondence be- 
tween the solution coefficient of expansion p, and the scalar potential ex- 
ponent A. So in effect, once one ffxes the three ffeld equation parameters 
{n, Vo,A}, there remains a two parameter family of these background so- 
lutions. Note also that these solutions are identical to those appearing in 
[14j . with the simple addition of the condition (|8]). For convenience, we shall 
denote by ''{g, (j), A){to,p,co,K}" a particular choice of the background solution 
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to the Einstein-Maxwell-scalar{„ vb,A} field theory; in using this notation, we 
presume that the conditions (P)- (niII) hold. 

One of the key properties of any of the background spacetimes (T" x 
{to,p,co,K} corresponding to the solutions ©-([H]) is the accelerated ex- 
pansion they exhibit, and the somewhat peculiar causal structure which con- 
sequently characterizes them. In particular, one finds that if one fixes a time 
to and the corresponding Cauchy slice T"^ in a background spacetime with 
expansion parameter p, and if for any point q G T"^ one considers a pair of 
coordinate balls i?e-«£o(5) ^^^1 B^e-i^ioiq) in T"^ for the characteristic length 
£o ■= then the causal future of i?e-K^o(^) contained in the future do- 
main of dependence of i?3e-K4(g); in terms of standard notation (see, e.g., 
Wald [IS]), one has 



The basis for this result is the fact that, if one considers any future causal 
path with starting point (g,to) on the to Cauchy surface T" , and if one 



the path can stray from q on T"^, a straightforward calculatioiu (see also |14j ) 
shows that it is bounded from above by Eq. Hence no causal path starting 
inside Be-i^g^lq) can reach a spacetime point for which there are inextendible 
past directed paths which avoid -Bse-'^fol^)' result follows. 

Relying on this result, we can spatially localize the study of the GFC- 
stability of our background spacetimes, since in analyzing the future causal 
behavior of the spacetime evolved from perturbed data in i?e-«^o(9')' ^® 
need not consider the influence of the development of any data outside of 
Bse-i^eoiq)- Note that there is a small simplification of the proof of our re- 
sults below if we work with an exterior ball of radius 4e~'^£o rather than 
3e~'^£o. Also, it is convenient to rescale the spatial metric in our background 
solutions by choosing the constant n = := ln[4£(to)]- Doing this, we have 
the slightly simpler causal condition J~^[Bi{q) x {to}] C D^[Bi{q) x {to}], 
which we can exploit in stating and proving our results here. 

^It follows from the expression (|6]) for the metric that any causal path 7(5) satisfies 
the condition -(7*)^ + e^'*(^)^P%7*7^ < 0, which can be rewritten as e^''(5y7'7-' < 
(^)^P(7*)^. One then calculates the projected displacement as J^^[e'^'^Sijj'^^^]^^'^ds < 



J+[5e-«,(t„)(g) X {to}] C D+[fi3e-^(to)(?) >< {^o}]. 



(11) 
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1.3 Initial Value Formulation of the Field Equations 

The statement of our results, as well as the proof, rely strongly on a formula- 
tion of the field equations as an initial value problem. The standard 
n+1 ADM-type initial value formulation is as follows: The initial data con- 
sist of a choice of (i) a spatial manifold S", (ii) a Riemannian metric hab and 
a symmetric tensor Kcd on S" which together comprise the gravitational ini- 
tial data, (iii) a pair of scalar fields (f and vr on which provide initial data 
for the scalar field, and (iv) a two-form B and one-form E on S" which make 
up the electromagnetic initial data. The initial data set (S", h, K, if, vr, B, E) 
satisfies the Einstein- Maxwell- scalar constraint equations (consisting of cer- 
tain components of the field equations if the following holdlj 

R - KijK'^ + (tiKf = vr^ + V'lpV^ip + 2V[^) + {EjE^ + ]^BijB'^), (12) 

V^K.i - Vi{iiK) = ifV,^ + E,B^ i , (13) 

ViE' = 0. (14) 

Here V is the Levi-Civita connection of h, R is its scalar curvature, and the 
indices are raised and lowered using h. 

Note that if we choose one of the natural t = const. Cauchy surfaces (say, 
t = ti) of the background solution {g, (j), A){to,p,co,K,}^ then the initial data on 
this Cauchy surface ish = e'^'^^^yPSijdx^dx^ , K = -pe'^'^t^^ {^)'^P'^6ijdx''dx^ , 

= f Inti - ^, vr = f ^, 5 = 0, and E = 0. 

Given a choice of initial data satisfying the constraint equations, one seeks 
a globally hyperbolic development of the data, which is a set {M"'~^^, g,(j), A) 
such that the following hold true: (a) {M"'~^^,g) is a globally hyperbolic 
spacetime, with M""*"^ diffeomorphic to S" x R_,_; (b) {M"-~^^ , g, (p, A) satis- 
fies the Einstein-Maxwell-scalar field equations (c) there exists an 
embedding i : S"' — )■ M""^^ such that 2(E") is a Cauchy hypersurface for 
(M""*"^, (yf), with first and second fundamental forms h and K, with (poi = (p 
and Vex0 o 2 = vr (for e± the future-directed unit normal vector field on 
z(S")), and with B = i*F and E = i*F{e±, ■) for F = dA. 

With small modifications (to generalize from the vacuum Einstein equa- 
tions to the Einstein-Maxwell-scalar field equations), the well-known results 

•^Here and throughout the paper, Latin indices run from 1 to n (space only) while Greek 
indices run from to n (space plus time). 
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of Choquet-Bruhat [1] (see also and of Choquet-Bruhat and Geroch [5] 
guarantee that for any smooth set of initial data satisfying the constraint 
equations f|T2|) - f|T4|) there exists a globally hyperbolic development; more- 
over, for such data there exists a maximal globally hyperbolic development 
(MGH-development) unique up to isometry, which is maximal in the sense of 
containment (with appropriate isometry map). The existence and uniqueness 
of MGH-developments plays a crucial role in the statement of our results, and 
in the proof of GFC-stability. 

While it is useful to state our main theorem (below) in terms of initial 
data sets of the form (E", h, K, ip, n, B, E), in carrying out the proof of our 
results we are led to work with modified specifications of initial data sets, 
which include quantities such as g^i and 1 + (^oo- Inclusion of these quantities 
is closely tied with the need to control gauges in the analysis, as we see below. 

1.4 Main Results 

The standard idea of a stability theorem is that one fixes a solution of the 
field equations, noting certain important properties of the solution, one con- 
siders certain classes of perturbations of the solution, and one shows that 
the properties of interest remain true for the perturbed solutions. As a con- 
sequence of the localized character of the causal structure of the expanding 
solutions under study here (see Section [TT^ . following [H] we state our main 
theorems here in a slightly different way (which effectively leads to slightly 
stronger results). We consider sets of initial data (S", h, K, (p, tt, B, E) for an 
Einstein-Maxwell-scalar {„ vb, A} field theory which in local regions are small 
perturbations of local initial data for one of our background solutions, and 
proceed to prove that the future development of the data restricted to a some- 
what smaller region has the desired properties (causal geodesic completeness, 
etc.). With our results stated this way, they apply to solutions which may 
only locally be a small perturbation of one of the {g, 0, A){to,p,co,K} solutions, 
or may be a perturbation of one of them in one region, and a different one in 
another region. This allows the results to hold for solutions with unrestricted 
topologies (unlike the background solutions, which are assumed to have the 
topology T" X R+). 

To measure the degree to which initial data for the perturbed solutions 
locally deviate from that of the background solutions, we need to work with 
a set of norms. Since the proof here depends on control of these norms via 
energy functionals, we are led to work with Sobolev norms; since the anal- 
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ysis is essentially local, we work with local Sobolev norms. In particular, 
for an open set ?7 C diffeomorphic to a ball in R" and therefore cov- 
ered by Euclidean coordinates (x^, ...a;"), for a tensor field \l/ on Y7^ with 
coordinate-basis components and for a non-negative integer m, we 

work with Sobolev norms defined as follows 

/ \ 1/2 

Inn „ 

\h, - ,iq=l ji, - ,jr=l \a\<m -^^(^^ 

Here the collective multi-index notation "9"" is used for the partial deriva- 
tives, all of which are calculated using the coordinate basis. 

In comparing (locally) a given set of initial data (S", h, K, ip, vr, B, E) for 
a perturbed solution with the data of a background solution, it is useful 
to find the "closest" background solution for the comparison. Presuming 
that the parameters n, Vq, and A have been chosen — thereby fixing the field 
theory and also thereby fixing (via and ( ITU]) ) p and Cq — it remains to 
determine to and k. As discussed above, it is convenient to choose (as a 
scaling) kq := ln[4£(to)]- Hence, one needs only to determine to. 

The idea for determining the appropriate choice of to is based on equation 
([7]), which (for a given A and Cq) gives the time dependence of the background 
scalar field 0. Roughly speaking, to determine to one calculates from the 
given (perturbed) data a local average of the scalar field, and then setting 
4>{to) equal to this average and inverting ([7]), one obtains to. More precisely, 
one chooses an open set f/ C S", together with a diffeomorphism ( : U ^ 
-Bi(O) C R". Then, one calculates (ip) := ^ Ibi{o) ^° "^^ere a;„ is the 

volume of the unit ball in R" with respect to the Euclidean metric. Finally, 
one sets 

-(A((^) + co) 



to := exp 



(15) 



We note that for a given set of initial data, this procedure for mapping 
to a comparison background solution depends only on the choice of the open 
set U and on the choice of the map C, \U — )■ i?i(0); once that choice is made, 
to (and therefore the comparison solution) is uniquely determined, regardless 
of whether the initial data is indeed close to a background solution. To 
simplify the discussion below, we use the notation 9{(/^^}(E", /i, K, ip, vr, B, E) 
to denote the map taking the indicated data set to to, as defined above. 

We are now ready to state our main theorem: 
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Theorem 1. Let {'E'^,h,K,ip,7T,B,E) be a set of initial data satisfying the 
constraint equations n^) -P^ for a fixed choice of the Einstein- Maxwell- 
scalars^n,Vo,X} i^eW theory. There exists an e > (depending only on n and 
p ) such that if for some open set U <Z Y7^ and for some diffeomorphism 
C : f/ — )■ -Bi(O) C R" the data satisfy the smallness condition 

||e-2'^«/i - 5\\h^o+^{u) + We'^'^HoK - p5||/^-o(r/) 

2 

+ - (V')ll/f'"0 + l(C/) + PoTT - -^||h"'0(1/) 

+ ||-E'i||/i''"o(c/) + ||-SjA:||_H'™o(c/) < e, (16) 

i i,k 

with to = Q^jj^^iYT^ K.ip.T: , B , E) , with kq chosen as above, and with mo 
the smallest integer satisfying mo > n/2 + 1; then the MGH- development 
(M"+\c/,0,A) of {Y^,h,K,ip,T^,B,E) has the property that if i : TP ^ M 
labels the embedding corresponding to the initial data, then all causal geodesies 
starting in i{C^"'^[i?i/4(0)]} are future complete. 

This theorem shows that for a data set which in a local region is suffi- 
ciently close to data for one of the background solutions, geodesic complete- 
ness holds in the development of the data on a specified subset of that local 
region. If this holds for data in a neighborhood of every point in E", then 
clearly the entire MGH-development of the data set is future geodesically 
complete. It follows as a special case that if one chooses a Cauchy surface 
T"p for one of the background solutions {g.cj), A)i^to,p,co,K}, and if one consid- 
ers sufficiently small Einstein-Maxwell-scalar {„^yQ_A} held perturbations of the 
data on T"^, then the MGH-development of that data is future geodesically 

complete; hence the solutions {g^cp, A)i^to,p,CQ,K,} are all GFC-stable. 

One might ask if, in addition to the property of future geodesic com- 
pleteness, the MGH-development of a set of perturbed data has the property 
that its fields in some sense approach those of the corresponding background 
solution. This is in fact the case, in a certain weak sense: 

Theorem 2. Let {TP, h, K, ip, tt, B, E) be a set of initial data satisfying the 
constraint equations in^) -[l4\) for a fixed choice of the Einstein-Maxwell- 
scalar{n,Vo,M fi-^^d theory, and also satisfying the e-smallness condition ( 176]) 
from TheoremUl Let {M^^^,g,(j),A) denote the MGH-development of this 
data. There are constants t_ G {0,to), a > 0, and > for all non- 
negative integers m, there is a smooth map : (t_,oo) x B^/q{0) — )■ M^^^ 
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which is a dijfeomorphism onto its image and satisfies ^(to,?) = >■ ° (^^{q) 
for q G i?5/8(0), and there is a Riemannian metric Hah on Br^j^iQi) such that 
the following are true: 

All causal paths that start in HC~^[-Si/4(0)]} remain in Imagej^'} for all 
of the future. 

Letting \\ ■ H^m denote the C"* norm on B^i^{Q), letting {g, (p, A) denote the 
pullback of the MGH- development fields via and letting (g, 0, A) denote 
the corresponding background fields, we have, for t > to, the following decay 
estimates: 

U{t, •) - 4>mC^n + \\{tdtCt){t, •)) - {tdtm)\\crn < ^^mit/to)-" , (17) 

ll^.llc"^ = \\d,Ao-doA,\\cm<Kme''°{t/tor{t/to)-'-^, (18) 
WBijWcm = \\d^A, - d,A\\c^ < Krr,e^^^^{t/tofP (tAo)"'"", (19) 
\\(l + goo){t,-)\\crn + \\tdtgoo{t,-)\\crr^<i^m{t/to)-\ (20) 

~ (^_2^)p+ i ^^'^^'^'ll^'" + \\^^tij9oiit, •))l|c- 

<t^m{t/to)-% (21) 

\\{t/to)-''Pe-^'''>g^,{t,-)-H,,\\c^ 

+ \\{t/toy^Pe-^''^'itdtgij{t, ■)) - 2pH,,\\cr^ < «„(t/to)~" , (22) 

\\it/to)''Pe''^'g'\t, ■) - W^\c^ < Kmit/to)-" , (23) 

||(^/^o)-'^e-2""ti^^,(t,■) -pH,,\\c^ < K^{t/to)-r (24) 

Here '-/ju are the (lowered index) Christoffel symbols for the metric H on 
i?5/8(0), and Kij is the (evolving) second fundamental form for the hypersur- 
face ^5/8(0) X {t}. 

We remark that while most of these inequahties clearly indicate decay, 
two of them appear not to do so: the second and the third, involving electro- 
magnetic fields. We note, however, that the electromagnetic fields are vector 
components relative to coordinates in which the metric is expanding. If one 
considers locally measured fields (factoring out the expansion), then these 
fields do decay. 

We also remark that these decay results do not show that the develop- 
ments of the perturbed data sets decay to the original background metrics 
directly. Rather, one obtains decay only if one adds a diffeomorphism, and 
also adds the fiducial metric H on i?5/8- It may be that sharper and more 
direct decay rates can be proven. We do not pursue this question here. 
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Finally, we note that while the results we prove here are generalizations 
of Theorem 2 in [T3], we have not gone on to prove a generalization of 
Ringstrom's Theorem 3, which applies his Theorem 2 to prove global sta- 
bility for a class of locally spatially homogeneous spacetimes. Such results 
likely could be obtained for spatially homogeneous spacetimes containing 
electromagnetic fields; we do not, however, consider that issue here. 

1.5 Outline of the Proof 

A key feature of the proofs of Theorems [1] and [2] is the spatial localizability 
of the analysis, noted above. This allows the analysis to be carried out 
independently on each open set U satisfying the hypotheses of Theorem [TJ 
However in order to avoid working on regions with free boundaries, it is useful 
to patch the data set on U into a set of background solution data on T" \ U, 
and then study the development of the patched-together data on T". The 
expansion behavior of the background solutions as well as their perturbations 
guarantee that the development of the data on sufficiently small subsets of 
U is independent of the externally patched-in data; hence the proof can be 
done via analysis on these patched data sets. 

In general, the patched initial data sets violate the constraints in an 
annular region around U. It is thus necessary to formulate a global stability 
analysis that works for initial data sets which violate the constraint equations. 
We begin to set up such an analysis in Section [2J To start, we modify 
the field equations by introducing gauge source functions and Q, and 
using them to hyperbolize the field equations. The gauge source function C 
replaces the contracted Christoffel symbols of the perturbed unknown metric 
by that of the background metric we perturb around and therefore is related 
to the "wave coordinate gauge", while Q, if it is zero, corresponds to the 
Lorentz gauge of electromagnetism. For handling a global stability problem, 
hyperbolization of the field equations is not enough. We also need to add 
to the field equations some correction quantities that are expressed in terms 
of the gauge source functions. The purpose of these correction terms is to 
partially decouple the field equations to linear order in the field perturbed, 
and to insert damping terms in these equations. These features are helpful 
in proving stability for the modified system. Since we are interested in global 
stability of the original field equations, we reformulate the equations in such a 
way that the gauge source functions satisfy a system of hyperbolic equations 
which admit zero as a solution, and so that it is possible to prepare initial 



13 



data for the modified equations so that the gauge source functions and their 
first order time derivatives vanish initially. Thus, the gauge source functions 
vanish identically, and by proving global stability for the modified equations, 
we obtain global stability for our original field equations. 

We express our reformulated equations as PDEs for variables which are 
essentially the differences between the perturbed fields and the corresponding 
background fields. More specifically, we work with u := 1 + goo, Ui := goi, 

:= (f) — (j), Ao, and Ai, all of which vanish for the background solution, 
and also 'jij := (t/to) "^^gij which is time independent for gij = gij. We 
further make a change of time coordinate from t to r, such that tdt = dr, 
by defining r = In (t/to)- The reason for doing this is to eliminate the time 
dependences of the background fields. Thus we obtain a system of equations 
f l65|l - fITOl) such that each equation is in the form of a hyperbolic equation 
with dissipation and dispersion, plus some extra terms. Note that in the 
modified system, the equations for u, ip and Ao are decoupled to first order 
from the equations for the other field variables. 

In Section [31 we define a sequence of energy functional for the field 
variables, and we specify the bootstrap assumptions which we use to prove 
global existence in r. The bootstrap assumptions state that for all r G / = 
[0, s) for some unspecified s, we have solutions to the field equations, and 
the energy functionals are controlled by some small number e. Note that 
the specific bootstrap assumptions we use here are not optimal; they are 
chosen because they are sufficient to carry out the global existence argument, 
and because they make the estimation of the nonlinear terms (following the 
algorithm for estimates introduced in |13j) applicable in our case. 

One of the key tools for proving global existence is the set of differential 
inequalities f ll40p - fll44p for the energy functionals. We derive these in 
Section [31 and we also show that as a consequence of the hierachical structure 
of the equations, the differential inequalities exhibit an hierachical structure 
as well. Relying on these inequalities, we are able to show (in Section [4]) that 
the bootstrap assumptions can be improved. That is, combining the the 
bootstrap assumptions with f ll40p - fll44p . we can show that a more rigorous 
version of the bootstrap assumptions holds. Bootstrap improvement is then 
used (also in Section [4]) to verify the "open" portion of an "open-closed" 
type argument which show that indeed, the interval / on which solutions 
exist extends to s = oo. We thereby prove global existence. 

The remaining work in proving our two theorems is first to show that the 
geodesies which start from a subset of the domain U x {to} are complete (we 
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do this in Section [S]), and then to verify the asymptotic expansions (i.e., the 
decay results) f ll89p - (11961) for the fields in the MGH-development of our 
initial data, as stated in Theorem [21 we do this in Section [61 We conclude 
our proof of the main theorems in Section [71 and make concluding remarks 
in Section [81 

2 Field Equation Reformulation 

As noted in Section [H as a consequence of the accelerated expansion of the 
background solutions, we can carry out the analysis in spatially local regions, 
which for convenience are each patched into a set of data for a background 
solution on T" (generally different for each local region). We defer discussion 
of the details of the patching to Section [71 For now, we presume that the 
patching has been done, and that we are consequently working with a set 
of data on T" which is (in an appropriate sense) a small perturbation of a 
set of background field data, but does not necessarily satisfy the constraints 
everywhere. The MGH-development of the patched data on T" consists of 
fields defined on the spacetime manifold T*^ x I for some interval /. Working 
on T" X J, and working with fields which are small perturbations of the 
background fields, we can always choose coordinates (x*,t) with global 
periodic spatial coordinates. 

The aim of the reformulation of the field equations ([ID-(|1|) we carry out 
here is to replace them by alternative equations which are manifestly hyper- 
bohc and lead to well-posedness even if the constraints are not satisfied, and 
also to make sure that the reformulated equations can be used to develop 
a set of energy functionals which are controlled in time and lead to global 
existence. The replacement equations are obtained by adding terms to equa- 
tions which can be made to vanish via gauge choice. We carry out the 
reformulation in two steps. In the first step, we add gauge terms which re- 
sult in equations which are manifestly hyperbolic for the fields g^j_ui (f) and 
(with components defined using the coordinate basis {x\t) defined above). 
In the second step, we rewrite the field equations in terms of the variables 
u, Ui, 'jij, ip, Aq, and Ai in a certain semi- decoupled form which is very useful 
for the analysis. 
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2.1 Reformulation I: Hyperbolization 

The gauge functions we use to carry out the first stage of reformulation of 
the field equations (II])- (H]) are 

:= f'^ - r'^, (25) 
Q := V^A„ (26) 

where := \g^^g^''{dagi3v + dpg^a — d^ga/s) is the contracted Christoffel 
symbol for the metric gai3, where := T^{g) = is the same for the 
background metric g, and where V is the covariant derivative compatible 
with gai3. Note that is not covariant. However, we define V^, := g^^V^, 
and also Vf^V^ := df,V^ - T^^V^. 

Using these gauge quantities, we define modified versions of the Ricci and 
Faraday tensors 

R^iu = Rfiu + ^{pViy), (27) 
F^, = F^, + g^,{g-V.,A''), (28) 

and we also define certain arrays of correction terms 

V, 





2p 




■~ T 




:= -9 




■= -9 



(29) 
(30) 



V^F^. + ^goA'DyA^ -G). (31) 

We then construct the following set of "gauge-modified" Einstein-Maxwell- 
scalar{„ A} field equations 

R,u + M^^l = d,ct>d,ct> + -^V{<P)g,, 

+ Ffj_(jFu — -^j— —g^jjFpijF'^ , (32) 

g"'^dc,dp(f) - F^(9^0 - V'{(j)) + M^*'^ = 0, (33) 
V^F^, + M^^l = 0, (34) 

noting the following properties: (i) If the gauge quantities and Q vanish, 
then this system fl32l) - fl34l) is satisfied if and only if the Einstein-Maxwell- 
scalarin^Vg^A} system is satisfied, (ii) The gauge modified system ( l32|l - 

( l34l) is manifestly (second order) hyperbolic for the fields {g,(j),A). 
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In view of these two properties of the gauge-modified system, our goal 
now is to show that for any given set of initial data0 (h, K, Lp,7i, B, E) (on 
the torus T") for the Einstein-Maxwell-scalar{„^vb,A} system there is 

a corresponding set of initial data {g, dtg, (p, dt^p, A, dtA) for f l32|) - flM|) which 
enforces the condition that the gauge quantities and Q vanish for as long 
as the solution exists. We show in the next subsection that this can be done. 

2.2 Initial Data and Gauge Choice 

There are many choices of the data {g, dtg, (f, dt(p. A, dtA) which are consis- 
tent with a given set of initial data {h, K, (p, vf , B, E) specified on T"; this 
is a manifestation of the gauge freedom in the Einstein-Maxwell-scalar field 
system. The following choice restricts this freedom, and in so doing, it en- 
forces the condition that the gauge functions and Q vanish at the initial 
tim^ to: 



9tjito, •) = Hdi, dj), gooito, ■) = -1, 9ojito, 


■) = o, 


(35) 




dtgijito, ■) = 2K{di,dj), 


(36) 




dtgoo{to,-) = 2t%to,-)- 


2tr^, 


(37) 


dtQoiito, ■) = 


-huT' + h''{2d,h^i-dikj) 


(to,-). 


(38) 




(f){to, ■) = (p, dt(t){to. 


■) = ^, 


(39) 




Ao(to,-) = 0, Ai{tQ,-) = 


Ai{x), 


(40) 


dtAo{to, ■] 


= h'^d,A,{to,-), dtAito,-) 


= —Ei, 


(41) 



where Ak{x) is a solution of diAj{x) — djAi{x) = Bij{x)^ 

The vanishing of the first time derivatives of and Q on (portions of) 
the initial surface follows not from further restrictions on the gauge choice, 
but rather from the constraint equations ( fT2|) -(fT^. The constraints are not 
satisfied everywhere (recall the consequences of the patching of the initial 
data) but they do hold on a subset S of the initial hypersurface. To see that 

^We use the tildes here to denote the fact that the data on T" has been constructed 
via patching of local data on U C S" to background data on T" \ U. 

^Note that T>o — corresponds to (l37l) . 2?; = corresponds to ((38)) . and Q = 
corresponds to the first part of (HTI) . 

^We choose an appropriate Ai{x) in Section [T] 
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dtV^^ and dtQ vanish on we first calculate the quantity G^y — T^i^, assuming 
that the gauge-modified field equations hold, rather than the original system 
(P-dl]); we obtain 

G,. - T,, = -V(^I?.) + \{V'^V,)g,, - M^^l + ^(^^"''M^fl)^?^,. (42) 

If we now let be the unit normal to the initial surface and let X be any 
vector tangent to the surface, and if we contract fH2|) with e_L and X while 
assuming that the constraints ffT^ - ffT^ hold on 5*, then we obtain 

- ]^{e^rX''{d^V, + d,V^) = 0. (43) 

Setting t = to, and noting that since Vi(tQ, ■) = we must have X'^d^V^itQ, ■) = 
0, it follows that 

dtV,{to,-) = OonSCT\ (44) 

Arguing similarly, but now contracting (H2l) twice with e± while assuming 
that the constraints hold at to, we obtain 

dtVoito, ■) = on 5 C T". (45) 

Finally, if we contract (IMI) with e±, then (setting t = to) we may use the 
constraint (fT^ to argue that e'^'V'^F^^ito, ■) = 0; combining this with the 
vanishing of V^(to,-), Q(to,-) and dtV^ito,-) (as shown above), we are left 
with 

dtGito, •) = on 5 C T". (46) 

Now that we have determined that for any choice of initial data we may 
choose a gauge so that Q, V^, and their first time derivatives vanish on 
the subset 5* of the initial slice, if we wish to show that Q and vanish 
in a spacetime neighborhood of S (within its domain of dependence) it is 
sufficient to show that these quantities satisfy homogeneous wave equations. 
To do this, we first take the divergence of equation fB21) . noting that 

V^T^, = -M[^l9,0 - '^(M]:^] +d^ig- V^A^)). (47) 

It follows that, presuming the gauge- modified equations have smooth solu- 
tions, we obtain 

g'^^d^dfsV^ + ""^d^Vp + ^d,g + 'V, \W^Q = ^ (48) 
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for smooth functions Q^"^ , , V^^ and W^. Similarly, taking the diver- 
gence of (El]), we obtain 

g'^^d^dpG + md^G + r^dj:)p + jg + l^v^ = o. (49) 

for smooth functions H'^, /"^, J and L^. This pair of equations together 
constitutes the desired homogeneous hyperbolic system for Q and V^. 

2.3 Reformulation II: Perturbation Variables and First 
Order Semi-Decoupling 

The system of PDEs fl52]) - flM|) for the field variables {g^u, ^fi, (p) together with 
the system fH5]) - fH^ for the gauge quantities {V^,Q) constitute a coupled 
PDE system which could be used to argue local existence of solutions. To 
be able to show long time existence, however, it is advantageous to replace 
the variables {g^u, A^, 0) by others which are closely related to perturbations 
of the background solution; following [H], we choose to work with u := 

1 + goo, Ui ■= 9oh lij ■= {t/to)~^^gij, ip ■= (l>-<P, Aq, and Ai. 

In deriving (from fl32|) -f lM|) ) the evolution PDEs for these new variables, 
we wish to segregate those terms which are linear in perturbations of the 
background solution (i.e., linear in u,Ui,'jij,'ilj, Aq, and Ai) from those which 
are higher order. Doing this, we obtain the following 



o 2Xp{np- 1) 
—dt^ ^5 ^ + ^00 



g^ d^d^Ui + —dtUi + 



t ^ ' 

np 2p 

-da., + - 

-g^^'d.d^A, + + ^Q^Ao + A',, 



-g^^d^d^^ij + —dtiij + {X{np - 1)^ - u) + A^ 



0, 


(50) 


Ui 




0, 


(51) 


0, 


(52) 


0, 


(53) 


0, 


(54) 


0. 


(55) 
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Here, the quantities A^^, A^,, A^ and AJ,^ are all of quadratic order or 
higher in the quantities that vanish for the background solution (such as 
u,Ui,il), Aq, Ai, their first derivatives, and the first derivatives of '^ijja- More 
specifically, one can write 



A' 



00 



A' 



Oi 



Aoo-2 
A0.-2 



FoaFi °" — 



1 



2(n- 1 
1 



-dmFaiiF"'^ 



2{n- 1) 



goiFasF'^^ 



A^ = A 

a: 



-2p r 



p. p.^ 

2(n 



1) 



dijFafsF"''^ 



A' 



a: 



bA 



A 
2 

2 



Om 



^gndpA, - jdtAo - ju{t>^A^ - gf^d^A,] 



P 



^go^{T^A,~g^-^,A^ 



(56) 
(57) 

(58) 

(59) 
(60) 

(61) 



where A^j, and A^ are written out explicitly in the work of Ringstrom; see 
equations (51), (52), (55) and (56) in [H], together with equations (82)-(87) 
and (92)-(93) in pjj. 

We note that the parameters n,p, and A appear in (I50l)-(l55l) because the 
new variables (u, Uj, 7jj, ■?/', Aq, Aj) are defined in terms of the background 
solution {g, 0, A){to,p,co,K}- We also note that the quantity gf^" as well as other 
metric and Christoffel quantities appearing in (jSOJ-dSSD may be viewed as 
functions of u, Ui, and ■jij. 

The calculations leading from the PDEs fl32|) -fp^ 
tially the same as those done in proving Lemma 3 in 
(^S^, we calculate as follows: 



to fl5Q|) - fl53|) are essen- 
To derive (l5l|) and 



= V'F^^ + M^^] 



g'^'idpF^ 



r7 p 



Tl^F^^) + d^{V>'A, 
2p 



- V^d^A^ - A^d^V^ + -^g,^{T^A^ - gP'^d.A, 

g'^Opd^A^ - g^'do.dpA, - T^{d^A^ - d^A^) 

- g'^'Kpd^A^ + g'^'Kpd.A, + d^igp^'d^A^) 



^We include the primes on these A' quantities here because below, we replace them by 
slightly changed A quantities without primes. 
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= gp'^dpd^A^ - t^d^A^ - A^dJ:>^ - gP'^K.d^A^ 



2p 



+ gP'^K^d^A, + {d^gndpA, + f ^oa(rM^ - g'^'dpA,) 
= gf'dpd^A^ - r'd^A^ A^dat^ - 2g''^V'i^d,A^ 
+ ^goo.{t^A,-gP"d,A„). 

For q; = 0, we obtain 

= V^F^o + 
= gP'^dpd^Ao - t^^d^A^ - A^^ot^^ - 2g'^''ri^^d,A^ 

+ ^g^{r'^A,-g^d,A,) 
= g'^'d.d^Ao - "fdtAo + - 2gP-T^^d,A^ 
- ^fM^ + + ^(1 + goo){r'A, - g<"'d,A,) 



g'"'dpd,Ao - "^dtAo 



t 



np{2p — 1) 



+ 



-2g'''^T^,^d^A^ + ^gP'^d.A, + ^«(f^A^ - g'^^d^A,) 



= g'^dpd^Ao 
+ 2 



n + 2)p 



np(2p-l) 
Ct^o ^0 



- gP'Tl^)d,A, + ^u{t^A, - g'^^d.A,) + ^-d^A^ 



r[A] 



For a — i^ we obtain 

= V^F^i + 
= g^'^dpd.A, - f'^a^A - A^.d.t^' - 2gf"'V%d,A^, 

+ ^^g^i{V^A^ - g^'^d.A,) 
= g'"'d,d,Ai - ^I^—^dtAi - ^diAo 



+ 2 



-jFoi + jgoiit'^A, - g^'^d.A,) - g'Tf^d^A, 
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The evolution PDEs fl50p - (133]) involve a number of factors of t~^. These 
can be conveniently removed by multiplying all of the equations f l50p - (|3^ 
by t^, by replacing t by r := ln(t/to) (so that consequently dr = tdt), and 
by setting tq := In(to); we then obtain the following system of PDEs, to be 
solved for {u{x,t) = goo{x, e'^^^^^*°^) = goo{x, e'^^^°) , Ui{x, r) = etc.): 



DgU + [{n + 2)p - l]drU + 2p[n{p - 1) + l]u 




-jdrip - 2Xp{np - + Aqo 


= 0, 


(65) 


DgUi + {np — l)drUi + p{n — 2){2p — 


l)Ui 




-e^+^'>[jd,^ + 2pg'^Tum]+Ao, 


= 0, 


(66) 


Ogjij + (np - l)d-r'jij + 2p {\{np - 1)^ - u) 'jij + Aij 


= 0, 


(67) 


2 

Dgip + {np - l)drip + 2{np - l)^) - -m + 

A 


= 0, 


(68) 


UgA^ + [in + 2)p - l]drAQ + np{2p - 1)Aq + 


= 0, 


(69) 


UgAi + [{n - 2)p - l\drAi + 2pe^^°d^Ao + A^,, 


= 0. 


(70) 



Here we define the hyperbolic operator \3g via 

Ug ■= -g'^^dl - 2e"+"o^°*9,9i - e^^^+^'^g'^didj, 

and we define the nonlinear remainder terms A^j,, A^, Ag and A;,_j via 

Aoo:=(l+^?°V« + e'(^+^»)A[,o, (71) 

Ao. := (1 + g'')dru, + e2(^+^»)A[,„ (72) 

A,, :=(l + ^7°°)9.7., + e^(^+-'')A^, (73) 

A^:=(l + ^7°°)5.^ + e2(-+-")A;, (74) 

Ae:=(l + (7°°)9.Ao + e2(-+^°)Al, (75) 

A,, := (1 + g'')drA, + e^^^^^") A^,. (76) 

Our stability analysis in the rest of the paper focuses on the evolution 
PDEs (!65l) - (!70l) . We note that, if we ignore the A terms in these equations 
and also ignore (for the moment) the dependence of the wave operator \3g on 
the metric, then we have the following semi-decoupled setup: (i) Equation 
( 1^ involves Aq alone, (ii) Equation (170]) involves only Ai and Aq. (iii) 
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Equations and (^E\i together form a coupled system for u and ip, in- 
dependent of the other variables, (iv) Equation involves only and 
u and ^z^. (v) Equation fl66|) involves and u and and 7ij, but not the 
electromagnetic variables. This semi-decoupled structure plays an important 
role in the analysis we carry out below. 

3 Energy Functionals, Bootstrap Assumptions, 
Estimates, and Differential Inequalities 

The key tool for proving global existence for solutions to a Cauchy problem 
for a hyperbolic PDE system is the set of energy functionals for the system. 
For a general (nonlinear) system, these functionals are neither canonically 
determined nor unique. However they may be obtained, the necessary prop- 
erties are (i) that (perhaps assuming certain a priori conditions on the fields) 
their future evolution is bounded, and (ii) that they control appropriate 
norms of the field variables. In this section, we obtain energy functionals 
for the PDE system (!65|) - (170|) . we state our bootstrap assumptions and de- 
rive estimates to control nonlinear terms (such as those contained in the A 
terms) and finally we derive the differential inequalities which, so long as the 
bootstrap assumptions hold, control the evolution of the energy functionals. 

3.1 The Energy Functionals and their Differential In- 
equalities for a Fixed Spacetime Metric 

The field equations (!65|) -( I70|) all involve the differential operator Dg, which of 
course involves the metric g. By definition of the field variables u,Ui, and 7jj, 
the metric g is closely tied to them, and consequently g evolves with them. 
This fact (a key feature of Einstein's gravitational field equations) must be 
taken into account in setting up the energy functionals and verifying their 
evolution properties. 

However, we start our discussion of the energy functionals by artificially 
decoupling the metric as it appears in Qg and elsewhere in coefficients of f l65|) - 
f lTU]) from the evolving fields u, Ui, and •jij. We do this, in this subsection, by 
fixing a (generally time dependent) spacetime metric g, basing Qg and the 
other coefficients on this fixed g, and treating u,Ui, and 7^- as independent. 
We recouple g and the field variables in the next subsection, with the help 
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of bootstrap assumptions. 

In defining the energy functionals, it is useful to treat the field variables in 
blocks, according to the semi- decoupling of the (linearized) evolution equa- 
tion, noted above. So we start by working with just u and ip, which we write 
collectively as the 2- vector 



u 



u 



The evolution PDEs fl65|) and fl68|) then take the form 

□„u + CdrVL + Ju + A = 0, 



(77) 



(78) 



where J and C are the constant matrices 
J = 



C 



2p[n{p - 1) + 1] -2Xp{np - 1) 
-2/A 2{np - 1) 

n + 2)p-l -8/A 
np — 1 

and A is the vector of nonlinear terms 

Aoo 

Since the only effect of the electromagnetic field on the evolution of u 
is via the A term, the form of the energy functionals we use for u and the 
calculation of their evolutions are formally very similar to what appears in 
Section 4 of Ringstrom in [T3]. Following that narrative, we first obtain a 
matrix T (see equation (76) in [TIJ) which diagonalizes the matrix J, and 



(79) 



then setting u = T 
C = T~^CT, we have 



u. 



T-^A, J = T~\JT = diag{\^,X+} and 



□„u + CdrU + Ju + A = 0. 



(80) 



We next define the base energy functional we shall use for u. Letting ccs, 
bi, and 62 be any set of positive definite constants, and using the notation 
gij ^ f,2{T+ro) gij ^ defin(§ 

£cs\^] ■= I: 



-g^^drii'dru + g'^diu'dju - 2ccs9 ii drii + hu' + h2V)dx, 

(81) 



^Here and below, the subscripts CS refer to the fact that u involves the metric "lapse 
function" as well as the scalar field. 
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where the superscript t on u* indicates the transpose of u. We then verify 
the following: 

Lemma 1. Let g be a fixed Lorentz metric on the spacetime T" x I for some 
interval I, and let u be a solution to equation ^8^) on T" x / for some choice 
of the constants p > 1 and n > 3 characterizing an Einstein-Maxwell-scalar 
field theory. There exist positive constants rjcs, Ccs, bi, 62 andccs (depending 
on n and p) such that if we define £cs using the last three of these constants 
( as in flHTj) ) and if we require 



then 
and 



+ <r/£5, (82) 
^£5[u] > Ccs I (5rfl*5rU + g'^diU^djU + u*a)dx (83) 

< -2vcs£cs + [ [-(5rii* + ccsu')^ + ^csMdx, (84) 



dr 

where Acs[ii] is a function quadratic in u and its derivatives, defined in 
equation (8^ below. 

Proof. Examining inequality ( l83l) , and using the choice of rjcs to ensure that 
(7°° is close to —1, we see that so long as c^^ < bi for i = 1,2, there must 
exist a constant (cs for which inequality (|83|) holds. To achieve c^^ < 6j, we 
set 

bi = \- + ccsCu, 62 = A+ + CCSC22, (85) 

where Cu and C22 are components of the matrix C defined above, and we 
note that so long as we choose ccs to be sufficiently small, this inequality 
follows. 

Now consider the time derivative of £^£5[u]; we calculate 
d£cs _ f f 1 ~t//^ , /^t\a - o ~t j~ o -t 



{-7:drU\C + C')drU-drU'Ju-drU'A 

- {{p - 1) + ccs)g''^diU^djn + c^sl^^up - ccsi^^Cdrn - C£5U*Ju 

- C£5U* A + biudrii + b2tjjdr'ijj + Acs [ii\}dx, (86) 

where 
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-d^-g'^ + ip-irg 



dig^'dru'drii - d^g'^drii'diU - ccsdrg^'^u'drU. 



(87) 



Note that the constant "(p— 1)" appearing in equation ( !86l) is labeled as H in 
the analogous equation in [H]; here, to avoid introducing an extra constant, 
we leave it in the form "(p — 1)". Using ( l85l) for bi and 62, and choosing ccs 
to be sufficiently small, we now obtain 



cs 



dr 



I + &)dr^ + ccs\dM^ - - 1 + ccsYg'^d^u'd.n 



Ccs 



Ccs 



}dx, 



+ 



{d^ir + ccsii)^ + Acs[n]]dx. 



It follows from Lemma 5 of [14j, that the matrix C + C* is positive definite. 
Thus, provided (7°° is close enough to —1 and provided ccs is chosen to be 
small enough, we obtain (IMj) for some constant rjcs- □ 



It is important to note that the differential inequality ( 18^ is useful for 
controlling the evolution of £cs only if we can establish estimates for A and 
'^cs- We do this below, using the bootstrap assumptions. 

Since we need to control derivatives of u as well as u itself, it is very useful 
to work with energy functionals which involve higher spatial derivatives. For 
that purpose, we define the following sequence of energy functionals (labeled 
by positive integers m) 



Ecs,m = Yl ^csld'^u] 

\a\<m 



(89) 



where 9" indicates the usual multi- index spatial derivative, of order \a\. It 
follows readily from this definition of Ecs,m and from Lemma [1] that the 
following differential inequality holds: 

Corollary 1. Presuming the hypotheses of LemmaUl Ecs,m satisfies 
dEcs,m 



dr 



< — '^VcsEcs,r. 
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|a|<m 



(90) 



Thus far, we have developed a sequence of energy functionals only for the 
pair of field variables u and ip (with artificially fixed g). To obtain a similar 
sequence of energy functionals for Ui,'yij,Ao, and Aj, again with the metric 
g fixed, it is useful to work with solutions of a model PDE 

DgV + adrV + f3v = F (91) 

for the scalar function v, where a > and /3 > are constants, g is a fixed 
Lorentz metric on the spacetime manifold T" x I (for n > 3), and F is a 
fixed function on the spacetime. We have the following: 

Lemma 2. Let v be a solution of equation f91\) on T" x I , with g, a, /3 and 
F as stated above. There are constants r^c, C > '^'^^ 7, > 0, depending on 
a and (3, such that if the given metric satisfies 

l^7°° + l|<^c, (92) 
and if we define the energy functional via 

£^[v\ ■=]- f (-g^O{drvY + g'^divdjv - 2jg'^%drV + 5v^)dx, (93) 

then S\;[v] bounds the following quadratic integral 

£v[v]>C [ [{drvY + g'^divdjV + L,3V^]dx, (94) 

where = if /3 = and = 1 otherwise, and £y[v] satisfies the differential 
inequality 

< -2r],Sv + / [{drV + -fv)F + Av[v]]dx, (95) 



where Av[f] is given by (Ol) below. If l3 = 0, then 5 = 7 = 0. 



Proof. There are two separate cases to consider, depending on whether /3 
vanishes or not in the model PDE ( |9T|) . If /3 > 0, then we see that if we 
choose 7 = a/2 and 6 = (3 + a^/2, and if we choose rjc so that g^^ is close to 
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— 1, then the inequahty holds. Calculating the time derivative of £v[v] 
we get 



dT 



{-{a - 'y){drv) + {6 - l3 - ^a)vdrV - /S^yv^ 
{{p - 1) + -f)g'^divdjV + {drV + jv)F + Av[v]}dx, 



(96) 

with 

AvM = - m~g'^)vd,v - 2^{d,f')vd^v - ^^-g""' d.vd^v - {dfg^'^drvf 

- {d,r)divd^v - ]^{d^g^^){d^vf + {^-d^g'^ + - 1)^*') d,vd,v 

- ^drg'\drV - 7(^?°° + l){drVf. (97) 

If we now substitute in the above formulas for 7 and 5, we have 

'^^^ = -\ f {a{drvf + aPv"" + (2(p - 1) + a)~g'' d,vdjv}dx 

(1/1 Zi J ryin 

+ / {{drV + -iv)F + Av[v]}dx. 

Combining this with the inequality f lM|) . we determine that for some constant 
r/c, the differential inequality fl^ holds. 

For the /3 = case, we choose 7 = and 5 = 0. We then readily check 
that, (for (7°° sufficiently close to —1), the energy inequality holds. As 
well, using the calculation (!96|) . we readily verify (for the /9 = case) the 
differential inequality (1951) . □ 

As with the energy functionals for u, discussed above, it is useful to pro- 
ceed from £y to a sequence of energy functionals involving higher derivatives 
of v: 

One then proves the following differential inequality result: 

Corollary 2. Presuming that the hypotheses of Lemma\^hold, the sequence 
of higher order energy functionals Ey^m satisfy the following inequalities: 

dEy m 



< -2r7e^v,m+ XI / {idrd''v+-fd''v){d''F+[Dg,d'']v)+Av[d''v]}dx. 

(99) 



\a\<m 
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Proof. If one differentiates equation f l9ip and then applies Lemma ^ the 
corollary immediately follows. □ 



To obtain the energy functionals for the fields Ui,'yij,AQ, and Aj, we 
now manipulate the evolution equations for these fields so that, with varying 
specifications of the function F and of the constants a and (3, these evolution 
equations (for each of the components of Ui,'yij, etc.) match the form of the 
model equation fjOT]) . For present purposes, we presume that 5^ is a fixed 
Lorentz metric, satisfying the condition \g^^ + 1| < rj, with rj a constant to 
be determined below. 

For Ui, we work with equation f l^ . If we set a = — 1 > and 
(3 = p{n — 2)(2p — 1) > 0, and if we set F equal to the negative of all except 
the first three termqj in (166!) , then we have an equation of the form (1911) for 
each of the components of Ui. It then follows from Lemma [2] and Corollary [2] 
that there exists a set of positive constants 75^, ^5^, 775^ and such that 
if we set 7 = 75^ and S = 5sh iii the expression ( 1^ for £y and if we then 
define 



(for m any non- negative integer) then the conclusions of Lemma |2l including 
the bounding condition ( 1M|) (with ^ = Csh) holds for -En^.O; and the differ- 
ential inequality (lUU]) (with appropriate choices of constants) holds for each 
Eu^^m- Furthermore, if we set 



then one verifies that similar sorts of bounds and differential inequalities 
(accounting for the summation) hold for these quantities. 

For 7jj, we work with equation ( E7|) . In this case, we set a = np—1, (3 = 0, 
and F = —Aij — 2p{X{np— — u)'yij; we then have an equation of the form 
(!9T|) for the components of 'jij. Lemma [2] now implies that for appropriate 
choices of constants 7x = 0, 6m = 0, rjM and (m, a bounding condition of 
the form (1941) and a differential inequality of the form (193]) hold for i^v[7ij]- 

^We note that, according to this construction, the function F inchides information — 
terms depending on 7^^, on Ai, etc. — which is not known. The formal derivation of the 
differential inequalities of the form (|84p or (|90|) still works, however. 




(100) 



|A|<m 




(101) 
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However, the most useful energy functionals for the metric field are 

EM,m ■■=\Y.11 i^v[d%,] + / e-^''^a>^{d%^fdx\ , (102) 

i,i \\\<m ^ ^ 

with oa = for |A| = and a\ = 1 for |A| > (and with a > a constant 
to be determined below, by flllSp ). Below in Section [3^ using bootstrap 
assumptions and consequent estimates, we obtain a differential inequality 
for each E_M.m which controls its evolution for all time. We also verify that 
the energies E_M^rn control Sobolev norms of ■jij and of its derivatives. The 
unfamiliar second (integral) term in the definition f ll02p is needed so that the 
energies EM,m do indeed control these quantities. We note that this term 
is also consistent with the condition that the sequence of energies EM,m all 
vanish for fields which are the same as the background fields [g, (j), A){to,p,co,n}- 
We proceed to the energy functionals for the electromagnetic fields. For 
the electromagnetic scalar potential Aq, we work with the evolution equation 
(EH])- It takes the desired form if we set a = (n + 2)p — 1,(3 = np{2p — 1), 
and F = — Ag. It then follows from Lemma |2] that for appropriate choices 
of constants ■ysp, 6sv, Vsv and (sv, a bounding condition of the form 
and a differential inequality of the form ( 195|) hold for 6^v[Ao]- The energy 
functionals for the scalar potential of primary interest are 

Esr,^ = ^v[d^M- (103) 

|A|<m 

We use bootstrap arguments below (in Section 13. 2p to obtain differential in- 
equalities which control Eg-p ^, and consequently control Aq and its deriva- 
tives. 

Finally for Aj, the electromagnetic vector potential, we work with equa- 
tion ( I70|) . To match the form of equation ( 191]) . we set a = {n — 2)p—l, (3 = 0, 
and F = —Ab^i — 2pe'^~^'^'^diAQ. It then follows from Lemma [2] that we have 
the constants jy-p = 0, 5vp = 0, ?7vp and Cvp! the main energy functionals 
for Aj are 

Evr,m = Yl Yl ^v[5'v4,]. (104) 

i |A|<m 

Corollary [2] implies that inequalities similar to fl99|) can be derived; control 
of Ey-p^rn is established in Section 13.21 Note that due to the vanishing of 7v-p 
and 6\;p, E^p^m does not control ||y4j||//m; since we do not need to estimate 



30 



||Aj||j:/m anywhere in this paper, we do not include in the definition of Exi-p^^n 
an extra term similar to the extra term in the definition of Ej^ ^- 

To establish appropriate differential inequalities for the energy function- 
al£^ Ecs,m, Esn,m, EM,m, Es-p^rn, and Ey-p^rn, and consequently to control 
their evolution and that of the fields u,ip,Uj,'yij, Aq, and Aj, we need to 
know something about the behavior of the function "F" corresponding to 
each of the fields. In the next subsection, we use bootstrap assumptions to 
establish this knowledge. Bootstrap assumptions also play a role in control- 
ling the nature of the evolving metric g which appears in the differential 
operator Dg as well as elsewhere in the evolution equations (!65|) -( [70|) . once 
we restore the relation between g and the evolving field variables u, Ui, and 
jij. We also discuss this in the next subsection. 

3.2 Bootstrap Assumptions and the Consequent Esti- 
mates 

The idea of bootstrap assumptions is that one assumes that the evolving fields 
satisfy certain conditions for t G / C R^, one uses those conditions to prove 
that certain estimates consequently hold (on /), and then one uses these 
estimates together with the evolution equations to argue that the solutions 
(satisfying the bootstrap assumptions) exist for all time t G R^. 

We use two bootstraps assumptions here. The first (which, following 
Ringstrom, we call the "primary bootstrap assumption") essentially says 
that the evolving metric stays Lorentzian for t E I. The second (called the 
"main bootstrap assumption" ) says that the energy functionals for the fields 
evolving from the perturbed initial data set remain small for t G /. We now 
state these more precisely, and discuss their immediate consequences. 

To state the primary bootstrap assumption, it is useful to work with the 
following notation: Let g he a real valued (n + 1) x (n + 1) matrix, with 
components g^^i, for /i G {i,0} and i G {l,...n}. We use g\, to denote the 
n X n submatrix with components gij, and we write (^b > if this matrix 
is positive definite. We use g~^ to denote the inverse of g and we use g^ 
to denote the n x n submatrix with components g'^^ , presuming that these 

""^"The subscript labels for these energy functionals Ecs.km, EsH.m, EM,m, Esv.m, and 
Evv.rn stand for, respectively, £apse and (Scalar field, iST^ift, A^etric, (Scalar T'otential, 
and Vector 7'otential. 
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inverses exist. As above, we let u[g\ := 1 + Qqq. To simplify notation herel"] 
rather than using Ui to denote g^i we write v[g] := {g^i, ■ ■ ■ , gon). As well we 
use u[g-^] := 1 + = 1 + {g^o - g'^u^Uj)-^ and v[g-^] := {g'^\ ■ ■ ■ , g^^) = 
(• • ■ , {g^^UiUj — goo)^^g''^ui, ■ ■ ■). We recall that g is defined to be a Lorentz 
metric if it is symmetric, has n positive eigenvalues and a single negative one. 
Noting that g is Lorentzian if u[g] < 1 and if g\) > 0, we call g a canonical 
Lorentz metric if these two inequalities hold. We use £„ to denote the set of 
canonical (n + 1) x (n + 1) Lorentz matrices. 

Now, given a symmetric positive definite nxn matrix M with components 
Mij, and given a vector w G R", we write (presuming the Einstein summation 
convention) 

\w\m:= {M,,w'w^y^\ 

If M is the identity matrix, we simply write \w\ := \w\id- We can now state 
the following: 

Definition 1. Let p > 1, a > 0, Ci > 1, rj & (0, 1), Tq and Hi be real numbers 
and let n > 3 be an integer. A function g : I x T" — )■ £„ satisfies the 
primary bootstrap assumption VBA{p^a,ci,ri,To,Ki,n} on an interval I if 

-\w\^ < e-^P'-^^'\w\l <ci\w\\ (105) 

\u[g]\<v, (106) 

for all w G R" and for all (r, x) G / x T", with kq = tq + Ki. 

We remark that in Section HI we determine to from the initial data and 
we set kq := ln[4£(to)] = In It then follows that ki = In 

It is useful to note the following immediate consequence of the primary 
bootstrap assumption (following Lemma 7 of [13]), which we use to control 
the inverse of the metric: 

Lemma 3. //(^ : / x T" — » satisfies the VBA{p^a,ci,r],To,Ki,n} on an interval 
I, then there exists a constant rjo G (0, 1/4) such that if i] < r]Q, thei^ 

\v[g-']\<2c,e~'^^-''''^\v[g]\, (108) 

^^The specific aim is to distinguish the vector m from the scalar u without using indices. 
i^We note that {v[g],v[g-^]) = goig°'- 
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Mg],v[g-'])\<'^c,e-'''^-^-Mg\\\ (109) 

\u[g-']\<Ar,, (110) 

^|^|2<,2p.+2.o|^|2^<^|^|2^ (111) 

hold for all u; G R" and all (r, a;) G / x T". 

To show that the energy functionals we have defined, so long as they 
are finite, control the fields and their derivatives, we need only the primary 
bootstrap assumptions. To argue this, we fix a time interval J, and assume 
that VBA{p^a,ci,ri,To,Ki,n} holds, with r] and a set as 

rj := min{r]o, r]cs/4:,'nsH/ 4:, Vm/ 4:, 'nsv/4:,r]);v/ 4:}, (112) 

a := ^ mm{p - 1, r]cs, VsH, Vm,Vsv, Vvv}, (113) 

where rjcs comes from Lemma [H while rjgy^, rjj^, rj^-p, and rjy-p are the corre- 
sponding constants which arise in the application of Lemma [2] to the other 
fields. We note (as a consequence of Lemmas [U and E]) that these constants 
a and t] depend only on the dimension n and the expansion parameter p. 
Invoking a convenient rescaling (which we justify immediately following the 
proof of Lemma [5]) in defining the following. 



Ere —P^arrn 7^ -2pr-2fi;„+2aT 77 

rp 4Ko+2ar 771 rp ^2aT rp rp —2pT—2KQ+2aT rp 

(114) 

we obtain a lemma which summarizes the desired (Sobelev) norm control of 
the fields: 

Lemma 4. Let the fields {u, m,, 7jj, ■?/', ^o. A} be specified on the spacetime 
manifold, I x T"' for some time interval I. If the metric g constructed from 
these fields satisfies the VBA{p, a, ci, rj, tq, ki, n} ( and therefore is Lorentzian) 
for rj and a defined by Ul^) and U13\) . respectively, and for real numbers 
Ci > l,To, and ni, then the inequalities 

e^'^miH^ + \\dr^P\\H-+e-^P-''>^-''^\\dMH^^^] < CE'JI^, (115) 
c^^MIh^ + \\dru\\H^+e-'^P-'^^-^^\\d,u\\H^^^] < CE'JI^, (116) 
e-^^^^'^^^iWu^y^ + \\drU,\\H^ + e~^P-''>^-''^\\d,u,\\H^n] < CE'J^^, (117) 
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g_2p.-2.o+a.p^^^^. _ 2pg,^U^ + e-^P-'^^-'^'Wdig.jWH^] < CE'J,'^, (118) 

^-2pr-2.o\\Qag^.y < c e'/,]^, < |«| < m, (119) 

e'^"[||Ao||H"^ + \\drAo\\H^ + e-'^P-'^^-^^\\d,Ao\\H^] < CE'J^^, (120) 
e-P^'^<^^''^[\\drA\\H^ + e^^P-'^^-'^^WdiAWnm] < CE'^^I^ , (121) 
/ioW on I, where the constants depend on c\, n and p. 

Proof. This result follows from the definitions of the energy functionals, and 
from the apphcation of Lemma [1] and Lemma [2l We note that in ap- 
plying these lemmas, we only use the conclusions regarding energies and 
norms (e.g., (l83ll ). not the conclusions regarding differential inequalities (e.g., 
(IM|) ). To reach these conclusions, it is necessary that the metric satisfy 
VBAi^a,ci,-n,To,K-i_,n}i but it is not necessary that the fields satisfy the field equa- 
tions (fi-'flZni)' □ 



To state the main bootstrap assumption, it is useful to define the following 
collective m^'^-order energy functional 

Em '■= Ecs.m + Esn,m + Ej^^^n + Esv,m + Eyp^rn- (122) 

We can then define the following: 

Definition 2. Let p > 1, ci > 1, tq, ki and e > be real numbers, let n > 3 
be an integer, and let mo > n/2 + 1 be an integer. Let rj and a be specified by 
and respectively. The fields {g, ip, A) satisfy the main bootstrap 
assumption M.BAip^a,ci,-q,To,Ki,n,e} on the interval I if 

1. ^ : / X T" ^ £„, : / X T" -> R an(i : / X T" -> R are C°°, 

2. g satisfies VBA{p,a,ci,v,To,Ki,n} on I, 

3. g, ip and A satisfy 

Ell^r) < 6, (123) 

for T E I. 

The primary use of the main bootstrap assumption is to control the (non- 
linear) A-terms which appear in the field equations f l65|) -f l70|) . Controlling 
these A terms, as well as certain commutator terms involving spacetime 
derivatives and the wave operator □, we can derive the differential inequali- 
ties (I140p - (I144p . which are crucial to our proof of global existence. We state 
the level of control of the A-terms in the following lemma: 
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Lemma 5. If the fields {g,il!,A) satisfy M.BA{p^a,ci,ri,To,Ki,n,e} on an interval 
I , then the following estimates 



(124) 
(125) 
(126) 
(127) 
(128) 
(129) 



hold on I, where the constant coefficients depend on n, p, m and Ci. In 
addition, the estimates 

||A£5[9"a]|Ui < Cte-^^Ecs,m , (130) 

\\Asn[d''uj]\\Li < Cee-''^Esn,m , (131) 

\\AM[d''h,j]\\Li < Cee-''^EM,m, (132) 

||A5p[5%]||li < Cee-''^Esv,m , (133) 

\\Ayr[d''A,]\\Li < Cee-''^Evv,m , (134) 



hold on I , for \a\ < m, where the constant coefficients depend here on n, p, 
m, Ci and on e~ 



-Kl 



Proof. In Section 9.1 of ^J, Ringstrom develops a systematic algorithm for 
estimating nonlinear terms similar to the A-terms discussed here, presuming 
that an appropriate bootstrap assumption holds. He applies it in Section 
6 of [14J to a set of quantities very similar to those of interest here. The 
proof of the present lemma involves essentially the same application; we 
omit the details. We do, however, note that the A^j, quantities here differ 
from the corresponding quantities in [H] only by terms which are quadratic 
in the electromagnetic fields. Since one determines that the electromagnetic 
fields have "good scaling properties" according to the classification scheme 
of the algorithm of [TB], the estimates calculated in |T3] for A^,^ hold here. 
Comparing the A^ quantity here and in [H], we find no differences, so the 
estimate from [H] applies. As for Ag and Ab^j, the estimates fll28l) - fll29p are 
again a straightforward application of the algorithm developed in [13]. The 
same is true of the estimates fll30p - (11341) . 

□ 
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We remark here that Ringstrom's estimation algorithm works optimally 
with the energy functional E^, with the scalings listed in f lll4p built into the 
definitions of Ecs,m., EsT-cm-, etc. This is the prime reason for working with 
these rescalings. 

The estimates just discussed, in Lemma |5l rely on the main bootstrap 
assumption, but do not depend on the fields satisfying any field equations. 
This next set of estimates do require that the field equations be satisfied, in 
addition to the main bootstrap assumption. 

Lemma 6. // the fields {g, ip, A) satisfy M.BA{p^a,ci,ri,To,Ki,n,e} on an inter- 
val I, and also satisfy the field equations ( 165|) -( [70|) on I , then the following 
estimates 



\\[n,,d-]n\\L^<Cee-'^^Ell' 
II [□„ 51^,1^2 < CeeP^+''°-2"^E;^/2 
||p5,51/^iilU^ <Cee2«o-2-EV2 

\\[Ug,d-]A4L^<Cee-'''^Eli^ 
I [□3,9"] A, 11^2 < Cee^'"+"°-2""E;^/2 



(135) 
(136) 
(137) 
(138) 
(139) 



hold on I, where < |q;| < m and where the constant coefficients depend on 
n, p, m, ci and on 6"*^^ 

Proof. Each of these commutators, when combined with the field equations 
f l65|) - fl70|) . produces a set of terms which includes one of the A-terms along 
with terms which are linear in the fields (or the metric inverses). The metric 
inverses are easily handled using Lemma [31 while control of the A-terms 
follows from Lemma E] above. □ 

With the estimates established in Lemmas [5] and [6] in hand, we may 
now derive the differential inequalities which control the growth of the en- 
ergy functionals, and consequently control the growth of the fields and their 
derivatives. 

Lemma 7. If the fields {g,i/j,A) satisfy MBA{p^a,ci,ri,To,Ki,n,e} o^s well as the 
field equations l[65\) - ([7^ on an interval I, then the following differential 
inequalities 

< -2aEcs,m + Cee-^'ETs,mE]l\ (140) 
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dEsH,m ^ „ p I /^pV2 (pl/2 , pl/2 \ , r^^^-ar pl/2 pl/2 
< -ZaJ^sn,m + '-^^5W,ml^£S,m + ^^,mJ + Eg^^^h^ , 

(141) 



(142) 

< -2aEsr,^ + Cee^^ E'JI^EK^ (143) 

< -2(p + 3a)Evp,^ + CEi/^^^E'Jl^ + C ee^ E^Jl^E^i^ (144) 

/ioW on I, where the constants depend on n, p, m, Ci and on e~'^^ . 

Proof. Fundamentally, these differential inequalities are all obtained by cal- 
culating the time derivative of the energy functionals, applying the field equa- 
tions, and then using the estimates of Lemmas [5] and [6] to control the non- 
linear terms. Since Corollary [T] and [2] (suitably applied) essentially contains 
the results of the time derivative calculations, the differential inequalities 
follow from that corollary combined with the appropriate estimates. 

We carry out the details here for Eg(^^ and E^^^, since the calculations 
for the others follow the same pattern, and differ little from those which are 
done in Section 7 of [H] 

Applying Corollary [2] to the field equation (1691) for Aq, and then (i) using 
the definition of Sobolev norms, (ii) using the estimates from Lemmas [5] and 
[6l and finally (iii) using the definition of E^^p^m, we obtain 



|«|<m 

+ I ^sv[d''Ao]dx 

|a|<m 

<-2r^svEsv,n^ + CY,i^\\^r^''M\2 + WAoh){WKh + m^^ 

\a\<m 

+ cY, iiA5PrAo]iii 

\a\<m 

< -2r]svEsv,m + C{\\drAo\\H- + ||Ao||//^)(|| Ae||//™ + Cee-2-E^/2) 
+ Cee~'''^Esp rn 
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< -2r)svEsv,m + Ce-^'^E^Jl^ee-^'^^El^^ + Cee-^^Esv,m. (145) 
If we then substitute in the relation Es'p,m '■= g'"^""^ Esv,mi we derive 



(146) 

the differential inequality f ll43p then follows. 

A similar string of calculations, starting from (ITOjl . produces ( I144p . Specif- 
ically, we calculate 



dE- 



vr,m 



I |a|<m 

+ 1 ^--^"^^ [-9"(2pe^+^°a,Ao + A,,) + [□„ a"] A] dx 

I |a|<m 

< —2riy'pEy'pm + CcC "''^Ey-pm 

+ 1 [-9"(2pe^+"°9,Ao + A,,,) + [□„ A] dx 

I |a|<m 

< —2riy'pEy'p„i + Cee ""^-Eyp.m 

+ \\drA,\\HACeeP^+^'-^''^EK^ + Ce^+^° ||9,Ao||h"^ + ||A,,,||^,™) 

< —2riy'pEy'p„i + Cee "''^Ey-p^rn 

+ Ce'^--''-^-"^-E'Ji,Jee-'^-El/' + E'J^J, (147) 
and then proceed to derive 

g2p.+2.o-2a.[2(p _ a)Evv,m + ^^^] < -2r^ypEvv,m + Cee-^^Evv,m 

I ^^„2pT+2Ko-3ar pl/2 pl/2 , ^ 2pr+2/to-2aT pl/2 pl/2 

from which (I144p follows. □ 

4 Global Existence 

A key step in the proof of our main theorem is showing that if a set of initial 
data satisfies appropriate smallness conditions, then the development of that 
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data exists for all future time. More specifically, as discussed in Section ll.4[ 
presuming that we have fixed a choice of the Einstein-Maxwell-scalar {„ vb, a} 
field theory, we consider an initial data set (S", h, K, vr, B, E). If this data 
set, in a region U E T,n, satisfies the smallness condition f lT6|) (which we 
justify below), then we calculate (y?) and to = Q{u,c}(X'"'yh,K,ip,Tr,B,E), 
thereby fixing a background solution (^, 0, y4){jQ p co.k} and a choice of time 
to- The choice of time fixes Kq := ln[4£(to)]. We proceed to smoothly glue 
{T,"',h, K,ip,'7i, B, E)\u into a copy of the background solution on T" \ U, 
producing data {h, K, cp, tt, B, E) on which agrees with the original data 
on U, satisfies the constraints everywhere except in a designated annulus, 
and satisfies the smallness condition ( !T6|) everywhere on T". We now show 
that the development of this data exists for all future time. 

The work of Section [2] shows that to prove global existence for the data 
(S", h, K, (f, 71, B, E)\u (or the data {h, K, (p, vf, B, E)) evolved via equations 
(IIJ-Q, it is sufficient to prove global existence for corresponding data, in the 
form {u,Ui,'yij,ip, Aq, Ai}, evolved via equations f l65|) - flTOj) . To carry this 
out, it is useful to first clarify the data correspondence: 

Definition 3. For a fixed choice of {n,Vo,X} — and therefore, correspond- 
ingly, a fixed choice of p and cq (see ([9]) and f lTOj) ) — let {h, K,ip,n, B, E) be 
a set of initial data on T", with corresponding values of (ip), to, and kq. The 
initial data for ([5^ - ( [7^ associated with {h, K, (p, vf , E, B) at t = is 

given by 

m(0,-) = 0, {dru){0,-) = 2np-2to{trk), (148) 
UiiO, ■) = 0, idrUi){0, ■) = Uoh'''i2dkhii - dihki), (149) 
7i,(0, ■) = kj, {dr%j){0, ■) = 2toKi, - 2phij, (150) 

7/^(0, ■) = ^- (<^), {d^tPM ■) = to^ - (151) 

Ao{0, ■) = 0, 9,v4o(0, ■) = toh'^diAjix), (152) 
A{0, ■) = Ai{x), drAiiO, ■) = -toEi, (153) 

where all the indices are defined with respect to the standard coordinates on 
Our global existence theorem for this data is as follows: 
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Theorem 3. For a fixed choice of {n, Vq, A}, let rj and a be given by M12\) 
and Ul!^) respectively and let tuq > n/2 + 1 be an integer. Let smooth data 
{h, K, (f, n, E, B) be specified on T", with corresponding values of (ip), to, and 
kq, and with associated data specified by ^148\ ) - U53\) . Assume that there is 
a constant ci > 2 such that 

-\v\'' < e-2"°7y(0,x)t;V < -\v\^ (154) 

holds for any v G R" and for x G T". There exist constants eo > and 
Cb G (0, |) (depending only on n, mo, p and Ci) such that if 

E'JUO) < cte, (155) 
holds for some e < eo, then there is a global solution to ^6^) - [JW. Further- 



more, 

miir) < 6, (156) 



holds for all r > 0, as do the primary bootstrap conditions l{105\) - (TO^. 



Proof. We fix a choice of e (possibly to be adjusted later), and we define A 
to be the set of s G [0, oo) for which there exists a set of fields {g, 0, A) on 
the spacetime manifold [0, s) x T" such that i) the fields {u, Ui, jij, ip, Aq, Ai) 
corresponding to {g, 0, A) satisfy the field equations fl65|) - (!70|) on [0, s) x T", 
and agree with the initial data (I148P - ( 1153^ at r = 0; and (ii), these fields 
satisfy the main bootstrap assumption AiBA{p^a,ci,ri,To,K.i,n,e} on [0, s) x T". 
To prove the theorem, we need to show that A is non-empty, relatively open, 
and relatively closed ( clS cl subset of [0,oo))@ 

To show that A is non-empty, we first note that, by construction, the sys- 
tem (|65|) - fITOl) is hyperbolic; hence, local existence for solutions of the initial 
value problem follows0 In particular, there exists an interval {Tmin,Tmax) 
containing r = on which there exists a unique (smooth) solution to fIBS]) - 
(ITOjl . To verify that this solution satisfies the main bootstrap assumption for 
some nonvanishing neighborhood of r = 0, we observe that since the initial 
data satisfy the conditions ^(0, ■) = 0, Ui{0, ■) = 0, and (11541) . it follows that 



^■^It is clear from its definition that A is connected. 

^""^We may also argue local existence for the initial value problem for (pS]) - (170)) by 
invoking the equivalence of this system, up to gauge, to the Einstein-Maxwell-scalar field 
system (as argued in Section [2]), and then relying on standard well-posedness theorems for 
the Einstein-Maxwell-scalar field system. 
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the primary bootstrap conditions fllOSp - f ll07p hold on some such neighbor- 
hood. Since the initial data also satisfies eUo{0) < fe, we see that fll56p 
holds on some neighborhood of r = 0. We thus verify M.BA{p^a,ci,ri,To,Ki,n,e}, 
and consequently determine that A is non-empty. 

To show that A is relatively closed, it is sufficient to show that if we 
assume that T G [0, oo) is contained in A, then the solution on [0, T) extends 
past T, and MBA{p^a,ci,ri,To,Ki,n,e} holds at r = It follows from the main 
bootstrap assumption J^BA{p^a,ci,r],To,K,i,n,€} and the field equations (|65l) - (ITOil 
that the fields {u,Ui,-yij,'ip, Aq, Ai) defined on [0,T) extend smoothly to T, 
with the metric maintaining its Lorentzian character. The well-posedness of 
the field equations with initial data at T (as argued above) then guarantees 
that the solution extends past T, and one readily verifies from the definition 
of MBA{p^a,ci,r],To,Ki,n,e} that if the solution extends and if the main bootstrap 
assumption holds on [0,T) then it also holds on [0,T]. 

It remains to show that A is relatively open. To do this, it is sufficient 
to show that, in the situation just discussed, the main bootstrap assumption 
holds not just on [0, T], but on some interval [0, T) with T > T. To show this, 
it is sufficient to show that an improved version of AiBA{p^a,ci,r],To,Ki,n,e} — i-e., 
a version with tighter estimates — holds on [0,T). 

The improvements we obtain on the 'PBA{p^a,ci,r],To,Ki,n} portion of the 
bootstrap assumption involve replacing the inequalities f ll05p - fll07p in Defi- 
nition [1] by strict inequalities. One does this by combining inequality (11230 
from the main bootstrap assumption with Lemma H] together with Sobolev 
embeddings to obtain control on the evolution of the metric components. 
The argument is identical to that in Theorem 4 of [14] . 

The key to improving the energy estimate (11230 is the set of differential 
inequalities in Lemma [7] which control the evolution of the energy functionals. 
Since AiBA{p^a,ci,ri,To,m,n,e} holds on the interval [0, T), these inequalities hold 
there as well. We consider first (I140p . Combining this with assumption (I123p . 
we have 

< -2aEcs,^, + Ce'e-E'Jl^^. (157) 
Integrating this first-order ordinary differential inequality, we obtain 

E'Js,m,ir) < e-^^E'Jl^^^iO) + \cre-'^^e' (158) 
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This proves that A is relatively closed, since it shows that if [0, T) C A, then [0, T] C A. 
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for all r G [0, T). Then using assumption (I123p to estimate Ej^^^^^{0), we get 

E'A"mo(^)<Ccs{c,e + e')e~^^/'. (159) 

Since the differential inequality fll43p for Esv,m{T) is essentially the same 
as that for Ecs,mo{'T)i we derive the same estimate for this quantity: 

EyimS'')<Csv{c,e + e')e-'^^/\ (160) 

To handle EM,mo{T), which satisfies fll42p . we first use the estimate fll59p 
to simplify f ll42p . and then use the integrating factor / = exp[—{e~"''^ — 1)] 
to solve the resulting first order ordinary differential inequality. We obtain, 
presuming that Cb < 1 (see Theorem 4 of 



EM,m,{r)<CM{che + t^)e. (161) 

To control Es-^^rnoi'T') , which satisfies fll4ip . we substitute the estimates 
derived above for Ecs,irio{j) and for EM,mo{.'T)i together with the bootstrap 
estimate assumed for Emoir), into the differential inequality (11410 . thereby 
obtaining 

— — — < -2aEsn,mo + Csn{Cb e + )^5w,mo- (162) 

We see that the right hand side of this inequahty is negative if EsH,mo exceeds 
a certain value; hence we obtain the upper bound 

Ei'Ljr)<^{cl'% + ^'^) (163) 

for this energy functional. 

The remaining energy functional to control is E\;-p ^moi'^)- Comparing the 
differential inequalities (I142p and (11440 . and noting from above (see estimates 
(11590 and (I160p ) the identical estimates for E^s^moi''') ^^"^ -^<sp,mo (''") 
see that the differential inequality for E'yp mo (''") is the same as that for 
-£'X,mo(T), except that it includes a decay term, — 2(p + 3a)i5vp,mo) and that 

~ ~ 1/2 ~ 

it does not include terms similar to Ce~"'^-E_v(,m and CE^g ^^Ej^ ^n- It im- 
mediately follows that 

Evv,m,{r)<Cvv{c,e + e^)e. (164) 
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We may now combine estimates ffT^ . ffMID . ffTMD . ffTB^ . and ffTM]) . to- 
gether with sufficiently small choices of and e (depending on the constants 
Ccs^ Cmi Csh, Csv and Cy-p) and conclude that 

^i/o'(^) < \^ (165) 

holds in [0,T). Thus we have an improvement on all estimates comprising 
M.BA[p^a,c]_,ri,To,Ki,n,e}- This proves global existence, thereby concluding the 
proof of this theorem. □ 

The global existence theorem just proven guarantees control of the energy 
functionals up to some chosen order uiq. In fact, we can show that the energy 
functionals to all orders m are finite: 

Theorem 4. For any solution of the field equations ( fg3]) - ( [70| j which satisfies 
the hypothesis of Theorem there exists a sequence of constants Cm such 
that for every positive integer m, the corresponding energy functional satisfies 

Ell\r) < Cm. (166) 

for all T > 0. 

Proof. It follows from Theorem [3] that the solutions under consideration sat- 
isfy the estimates of J^BA{p^a,ci,r],To,Ki,n,e} for all time r. Hence the differential 
inequalities f lUOj) - f Ulij) hold for all r > 0. 

It is useful to work with a set of three sequences of rescaled energy func- 
tionals; we set 

EsH,m '■= e "''^^'^EsH,m, Ecs,m '■= ^"''^^'^ Ecs,m, Esv,m '■= ^^^^"^ Esv,m- (167) 

Based on f ll40p . f ll4ip . and f ll43p . we find that the differential inequalities for 
these rescaled quantities take the form 

2aEsn.m + Ce--^"Eyi^{E]Jl^ + E'/,"m) 



dE, 



CS.m 



dr 

dEsv,m 

dr 
dEs^^m 



dr 



< 



< 



< - 
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If we now define the quantity 

£m '■= Ecs,m + Esn,m + EM,m + Esv,m + Ey-p^rn (168) 

tlien it follows from the above differential inequalities combined with (11421) 
and ( 1144P that £m{T) satisfies 

^ < Ce-^^/'£^ + CE'/I^Em,^. (169) 

We now set m = niQ. Noting the boundedness of EM,mo{'T) (see fll56p ) 
and noting the behavior of E^^ (r) (which we infer from the relationship 
(11671) between E]J^^^^{t) and -^^^^^^(r), and the relationship fll68p between 
^£5,mo(^) and £rno), we derive 

^ < Ce"-/^^„„. (170) 

It follows from this differential inequality that is boundeci^ and conse- 
quently that E]JI^^ < Ce-""^/^. 

^ Icy 

If we now substitute this estimate for E^^^^^ into f ll69p . we obtain (for 
arbitrary m) 

d£rn ^ Ce-^'^/^S^. (171) 
dr 

It follows that £m is bounded for all m, which implies that Ecs,m, E^M^^n, 
Esv,m and Ey-p^m are bounded as well. 

It remains to show that EsH,m is bounded. Substituting into (11410 the 
boundedness conditions just determined, we have 

dEg-j^ ~ ~ ~ 1 2 

— < —2aEs-H,m + Ce'"'^ Es-H^m + CEg^^. (172) 
ar 

For r large enough, the second term can be absorbed into the first term. It 
follows that EsH,k is bounded, since (11721) requires that it decay if it exceeds 
certain value. This proves the theorem. □ 



We note that the proof of Theorem [3] contains estimates of £m„ and of Ecs^m^ on a 
finite time interval Here, we show that the estimates hold for all r > 0. 
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5 Causal Geodesic Completeness 



With global existence for spacetime developments of certain types of initial 
data sets proven above, we argue in this section that causal geodesies with 
initial points contained in certain regions of those data sets are complete. We 
do this in two steps, with the first proposition establishing certain geometric 
properties of these developments. 

We note that since the behavior of causal geodesies depends exclusively 
on the metric, the discussion here is very similar to that in |14j . 

Proposition 1. Let {g,(j),A) denote a spacetime development (satisfying the 
field equations l[65\) - ([7^ ) of the sort constructed in Theorem Let •y be a 
future- directed causal curve contained in this development, with parametriza- 
tion domain [so,Smax), o-nd with 7''(so) = to ( to is given by ( [73]) ). If the 
e parametrizing the smallness of the data (as stated in Theorem\B() is suffi- 
ciently small (depending only on n, p and C\), then 7*^ > (for all s), and 
the projected spatial length of the path satisfies the condition 

/ [9^,ito.l^Wi'Y'^ds < d{e)i{to), (173) 

J so 

where d{e) is independent of •y, with d{e) — j- 1 as e — )■ 0, and where 7b(s) = 
(7-^(5), 7^ (s), ■ ■ ■ ,7"(s)). Furthermore, assuming thatj is future inextendible, 
we have 7°(s) — )■ 00 as s — )■ s^ax- 

Proof. Since 7 is a future-directed causal curve, we have the following in- 
equalities 

g^^uYY < 0, (174) 
goof + go^f<^. (175) 

It then follows from (11071) (one of the estimates comprising the primary 
bootstrap assumption, which we have shown in Theorem [3] holds for ((;, 0, A)) 
that for some constant 77 G (0, 1), 

|2(7o.7°7i < ^'/'|7T + r^'^\9o^if < v'^^H^'l^ + v'^^ci'e^P^^^^'''^''^6,,fy. 

(176) 

As a consequence of (I105P (another of the proven primary bootstrap esti- 
mates), the last term is bounded by i]^^'^gijYi'' ■ Combining (I106p and (I174p . 
we find that 

9i,yi' < c(r/)7V, (177) 
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where c{ri) — )■ 1 as 77 — )■ 0+. Again invoking f llOSp . we determine that 

< Cic(r7)e-2^--2''°7°7° = c,c{v){t/to)-'Pe-'^''j'^'. (178) 

Now, as a consequence of the energy estimate (11561) . the role of the energy 
functionals in controlhng function norms (see (IllSp ). and Sobolev embedding 
(recall the condition in Theorem |3] that mo > n/2 + 1), we have 

e-2'^°+"^||9,7.,||oo<Ce. (179) 

Consequently, we obtain 

\\{t/to)-'Pe-'^^gij{t, ■) - e-'^^'g,jito, OIU < Ca'^e, (180) 
where C depends only on n, p and ci. Combining this with f ll78p . we obtain 

le-'^^gijito, 7b)7¥ - {t/to)-'^e-'^'g.,Y¥\ < Ca-'ecic{r]){t/to)-'Pe-'^'iV ■ 

We note here that it follows from f lTT6|) . f lTTTll and f lT56D that|l!| r] in ffT06|l 
and (I107P can be replaced by Ce, where C only depends on n, p and Ci. 
Combining this observation with the two inequalities fll77p and (I18ip we 
produce the inequality 

e-^-'gijito, 7b)7¥ < t;'(e)(tAo)-'"e-2'^»7V, (182) 

where d{e) — )• 1 as e — )■ 0+. 

To show that 7° > 0, we consider fll75p . Combining fll07p and (I178p we 
derive 

\9oi¥\ < [e-'^^-'^'S'^9o^goJYV'^^'^'S^,f¥Y^' < e(e)|7°l, (183) 

where ^(e) — t- as e — )■ 0+. If we now make e sufficiently small (depending 
only on n, p and Ci), it follows from f ll83p and fll75p that 7*^ > 0. Combining 
this result with f ll82p . we obtain the estimate fll73p . 

To complete the proof, we assume that 7 is future inextendible and we 
suppose that 7° does not tend to 00. Since 7° > 0, it follows that 7° has to 
converge to a finite number. Examining f ll78p . we see that 7* must converge 
as well; hence 71, must converge to a point on T". But this contradicts the 
assumption that 7 is a future inextendible causal path. □ 
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The argument is as follows: From (|116p . (|117p and (I156p . we have |m| < < 
Ce-^^El^l^ < Cee-"^ and we have |ujf < \\uj\\%^ < C^e^P^+^'">-^''^ EsH^rn < 
(^2g2g2pT+2Ko-2aT_ jj: ^ jg gjj^g^}} enough, we can replace f] by e. 
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We proceed now to prove causal geodesic completeness: 

Proposition 2. For sufficiently small e (depending on n,p, and Ci), the 
spacetime development (satisfying the field equations l[U^) - (71^) of a set of 
initial data satisfying the hypothesis of Theorem\^is future causally geodesi- 
cally complete. 

Proof. We let 7 be a future-directed causal geodesic and we let {smin, Smax) 
denote the maximum range of its (proper) parameter. The geodesic equation 
for 7 takes the form 

7^ + r2.rr = 0; (184) 

for the purposes of our argument here, we work with the /3 = component 
of this equation: 

7° + err = 0. (185) 

We wish to estimate the second term of this equation. First, using the 
estimation algorithm Section 9.1 of [13] together with the energy estimate 
(11561) . we obtain 

|r°ol < C^e(p/t)e— , 
|^°,|<Ce(p/t)e^-+««-"^ 
|r°- - {plt)g,,\ < Ce(p/t)e2p-+2.o-a._ 

Consequently, for t large enough or e small enough, we have r^j-7*7-' > 0. 
Combining these estimates with (I178p . we conclude that 

|r°o7Vl + 2|r°,7°7i < Ce{p/t)e-^^\^y, (186) 

where C depends only on n, p and Ci. Based on these conclusions together 
with (11851) . we find that 

7° < Ce(p/t)e-'^rV = Ce(pA)(t/to)~rV, (187) 

for s > Si. We note here that 7° and t can be used interchangeably. 

Since 7° > (assuming e to be small enough), we may divide both sides 
of (11871) by 7° and then integrate, thereby producing 

In4f4 < Cep / t-\t/toy^°ds = Cep / t-Ut/t^y^dt < Cep/a, 

(188) 
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where we let Si be large enough such that 7°(si) > to if necessary. It follows 
that 7° is bounded away from 0. Hence we have 



s 



7°(5)-7°(so)= / i'{s)ds<C\s-So\ 



so 



for some constant C. Since 7°(s) — > oo as s — )■ Smax, we conclude that 
^max = c>o. Thus 7 is future complete. □ 



6 Asymptotic Expansions 

The results of Sections H] and O provide the bulk of the proof of our main 
result, Theorem [T] (we complete the proof below, in Section [7]). These re- 
sults say little about asymptotic behavior of the solutions of interest, beyond 
global existence and geodesic completeness. Here, we address the issue of 
asymptotic behavior, as a step toward the proof of Theorem [2l 

Proposition 3. For sufficiently small e (depending on n,p, ci, and mo), the 
spacetime development {g,(l),A) (satisfying the field equations / fg3]) - (7^) of 
a set of initial data satisfying the hypothesis of Theorem\^ has the following 
asymptotic behavior: There exists a smooth Riemannian metric H on T", 
and for every integer I > there exists a constant ai (depending only on n, 
p, C\ and I), such that for all t >to, the following estimates hold: 



^\nt+^\y + \\{tdt<P){t, 



jWc < o^iit/to) 



\\E.^\ci = \\d,A^ - d^A^Wci < aie^^^t/t^Y {t/t^)-^-^ 
\B^,,\\ci = \\d,A, - d.A^c^ < aie^^'{t/t^f^ {t/to)-''" 
\\{l + goo)it,-)\\c^ + \\itdtgoo)it,-)\\ci < aK^Ao)"" 



,1 / 
\j9oi{t, 



+ e 



(n -2)p + l 



H- 



jm 



Iji 



-2k() / 



Ic' + \\'tdt{jgoi){t,-] 
< aiit/to) 



\c- 



-2p 



~2ko 



{t/to 

{t/t^)-^P{tdtg,,){t,-)-2pH,, 
\\e^^-{t/tofPg^^{t,- 
\\e-^^^\tltor^nK,,{t,-) 



ijWC- 



c' < aiit/toY 
-W^\\c^<ai{t/t,y 
pHijWci < ai{t/toy 



(189) 

(190) 
(191) 
(192) 

(193) 

(194) 
(195) 
(196) 
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Here ^jim denote the (lowered index) Christoffel symbols associated to the 
metric H, Kij{t,-) are the components of the second fundamental form in- 
duced on the Cauchy surface {t} x T" by the spacetime metric g^^, o,nd \\ ■ \\ci 
denotes the C' norm on T". 

Remark. As noted above in Section 11.41 although it may appear from fll9UI) 
and fllQip that the electromagnetic field grows exponentially, this is essen- 
tially a coordinate effect, which reflects the use of coordinate bases relative to 
which the physical metric is expanding. If one examines the locally measured 
flelds (factoring out the metric expansion) then the electric fleld and mag- 
netic fleld decay as (t/to)"*-^'*'"''- In fact, one expects that a further analysis 
would show a stronger decay rate. 

Proof. The flrst four of these estimates, fll89p - fll92p . follow immediately 
from the energy bound in Theorem H] together with the norm bounds flllSp . 
(HinD, (HlOD and dm]). 

Based on the norm estimate (11 180 . on the relation := {t/to)~'^^gij = 
e'^^'^Qij, and on Sobolev embedding, we see that the Ci norm of e~'^'^° dr'^ij 
decays as e"""^. Thus there exist smooth functions Hij such that for every 
integer / > 0, 

||e-2'^«7i.-(r, ■) - H,,\\c, < aie'^^ (197) 

holds for some ai and for all r > 0. Combining this decay result with that 
for noted above, we obtain (11940 . 
To derive (I195p . we calculate 

It follows from (I USD that the Ci norm of the right hand side of the above 
decays as e"'*'^. Hence there exist smooth functions H^^ on T" such that 
(11950 holds for all r > 0. We also conclude from above that W^Hj^ = Si, 
which implies that Hij is a Riemannian metric on T*^ and H^^ is its inverse. 

To obtain the remaining estimates of this proposition, (I193P and (I196p . 
we refer to the derivation carried out in Section 10 of [H]. The corresponding 
analysis for the Einstein-Maxwell-scalar fleld theory requires that we handle 
an extra term: FQ^Fi'^ — -^^^^—^gQiFapF"'^ . We readily verify that this term 
decays sufficiently quickly for the argument of Section 10 of [11] to apply 
here. □ 
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7 Proof of the Main Theorem 



The theorems and propositions of Sections HJ [5], and [6] essentially prove our 
main results, Theorem [1] and Theorem [2l Here we discuss a few of the 
remaining details. We note that the steps needed here to complete these 
theorems are similar to those carried out in [13] and [11]; hence, we leave out 
some of the details. 

Step 1: Construction of a global-in-time patch, related to data on U. 
Theorem [3] of Section Improves global-in-time existence for smooth initial data 
{h, K, (f, TT, E, B) specified on and satisfying suitable smallness conditions. 
As noted in Section |U a key step in using Theorem |3] to prove our main 
result. Theorem [H is to construct this smooth data on T" from the given 
data {h,K,{p,7!-,E,B) specified on U G S" and satisfying the hypotheses 
of Theorem [TJ To do this, we define a cutoff function /c(x) G C^{Bi{0)) 
such that fc{x) = 1 for < 15/16 and < fc{x) < 1, and we define 
{h, K,<^,TT,E,B) on T" as follows: 

hij = fcihij o X-') + (1 - fc)e^'''Si„ (199) 

K,, = f,{K,, o X-') + (1 - e^'^"^,,-, (200) 

^0 

(p = f^ifO X'^ + (1 - fc){^) 

-^^[{fc{vox-'))-{^){^)], (201) 

^ = /c(vrox-i) + (l-/,)-^, (202) 

E = MEox-'), (203) 
B = MBox-'). (204) 

Here to and kq are defined in the statement of Theorem [T|, and the last term 
in (120 ip is included so as to ensure that to defined in Theorem [3] equals that 
defined in Theorem [H 

To go from the initial data set {h, K, ip, tt, E, B) just defined to the asso- 
ciated initial data for the evolution equations ( 165|) - ( 170|) . we rely on (11481) - 
(11531) . Included in this transition is passage from the magnetic field B{x) to 
the vector potential Ai{x). To show that this can always be done, with the 
needed norm control, we note the following^ 



^®This lemma is very similar to Lemma 2.6 in |15j . The main difference is that our 
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Lemma 8. Suppose that B{x) is a closed 2-form on T" such that suppB C 
-Bi(O) and Bij G _f/'™(^T") for m > n/2 + 1. There exists a 1-form A{x) on 
T" such that (i) B = dA; and (ii) for any given e > 0, there exists 5 > 
such that if 



< 5, (205) 



then 



Aj||/^m+l(X") < £• (206) 



Proof. Applying the Poincare Lemma in the form of Problem 15-2 from [10] 
to the closed two-form B{x) on the ball -Bi(O), we obtain a 1-form A on -Bi(O) 
such that B = dA and Ai = tx^ Bji{tx)dt on -Bi(O). Recalling that Bji 
vanishes at the boundary of the ball -Bi(O), we see that A does as well, so we 
may extend A continuously to all of T", with it vanishing on T" \ -Bi(O). 

To control the norm of the vector potential, we now carry out a gauge 
transformation: We replace Ahy A := A + d(j), with taken to be a solution 
of the Poisson equation A0 := — diAi. It follows from a straightforward 
calculation that the components Ai of the gauge-transformed vector potential 
satisfy a Poisson equation with the source term being the divergence of B. 
Applying standard elliptic theory (e.g.. Theorem 3 on page 249 of [12]) to 
this equation (on the closed manifold T"), we obtain estimates which lead 
to the conclusion of this lemma (see the proof of Lemma 2.6 in [13] for 
details of these estimates, done on a bounded domain rather than on a closed 
manifold) . □ 

The hypothesis of Theorem [3] includes the condition that Emo (0) be small. 
To show that this holds for the data (/i, K., (p^ vr, B) on T" constructed as 
above (from data satisfying the hypothesis of the Theorem [1]), we need to 
verify that for any given any 6 > 0, there exists e > such that if f lT6|) holds, 
then E}no{0) < The arguments leading to the estimates -^'^^^^^(O) < Ce, 

~ 1 /2 ~ 1 /2 

^SH,moi^) — ^M,mo(^) — ^^^^ similar to those in Section 

11 of [Ilj. Here we show that Egp^^{0) and Ey.^^^{0) must be arbitrarily 
small, presuming that we choose a small enough e. 



version is specifically carried out for T". 
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To show that Ejp,^^{0) can be made arbitrarily small, since Ao{0,-) = 

0, we only need to show that ||(9t-74o(0, ■)||h™o = \\toh^^diAj{x)\\H"^o can be 
made arbitrarily small. It is therefore sufficient to show that \\Ai{ 
can be made small enough. To estimate ^^^^^(0), we need to estimate 
\\e-''°drAi{0,-)\\H^o and ||e-''oAi(0, Olk^o+i- Combining fl2D^ . ffTB^ . and 
([T6|l . we conclude 

\\e-''"drAi{0,-)\\H-o <Ce. 

Hence it remains to estimate ||e~'*°ylj(0, ■)||j|/mo+i, and for this it is enough 
to estimate ||74j(x)||/^'mo+i- From Lemma [HI we conclude that ||Aj(a;)||j:/mo+i 
can be made arbitrarily small presuming that ||i?ij(a;)||//"io is small enough. 

~ 1 /2 ~ 1 /2 

Hence Eg.p ,^^{0) and Eyip ^^(0) can be arbitrarily small given that we choose 

'"1/2 

a small enough e. In conclusion, we have a small enough Enlo (0). 

The remaining condition from the hypothesis of Theorem [3] that we need 
to verify is the estimate (11541) . Relying on arguments similar to those of 
Section 11 of [H], we determine from (fT6|) that (11541) holds for some ci > 2. 
Thus, noting that Ki only depends on p and that mo only depends on n, we 
see that Theorem [3] applies here if we assume e to be small enough (depending 
on n, p). As a consequence, we obtain a solution {g', ip', A') to - (1701) on 
(t_,oo) X T" for some < t_ < to- K follows from Proposition [3] that we 
also have the asymptotic estimates (11891) - (I196p . 

We verify now that the development of data on certain subsets of -Bi(O) 
solve the Einstein-Maxwell-scalar field equations. Transforming the field vari- 
ables to [g', (j)', A') (in accordance with the discussion in Section [2]), we verify 
that these fields satisfy the modified field equations fl32]) . fl33]) and flMl) . Fur- 
thermore, since the constraints (fT2l) - (IT^ are satisfied on -Bi5/i6(0) and since 
the initial data have been constructed so that V^ito, ■) = dtV^ito, •) = and 
G(to, ■) = dtQ(tQ, •) = 0, it follows from standard hyperbolic PDE theory to- 
gether with arguments from Section [52] that on the domain of dependence of 
-Bi5/i6(0) X {to}, (which we label £)[i?i5/i6(0) x {to}]) the functions {g', 0', A'), 
satisfy the original Einstein-Scalar-Maxwell system of equations ([I]) - ([3]). 

Assuming e to be small enough, we conclude from Propositions [1] and [2] 
(following a line of argument similar to that found in Section 11 of [14j) that 
causal geodesies starting in i?i/4(0) x {to} are future complete, and that 

55/8(0) X (t_, OO) C Z}[529/32(0) X {to}], (207) 

where we increase t_ if necessary. Hence we have a global-in-time spacetime 
which contains the development of data on a subset of the set f/ G S" (from 
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the hypothesis of Theorem [T]) . Within this development, causal geodesies are 
complete. 

Step 2: Construction of local-in-time patches, away from U . To obtain 
a Cauchy development of the initial data on all of S", we also need patches 
that are developments of initial data on subsets of that are not contained 
in U . Such initial data does not generally satisfy any smallness condition. 

To work with data of this sort, it is useful to define a reference spacetime 
metric on x oo): We set 

g:={l-f,o x){-de + h) + if, o x){Id X xYg', (208) 

where h is the Riemannian metric on S" included in the initial data, g' is 
the solution metric from Step 1, /c is the cutoff function defined in Step 1, 
and X is a coordinate map. As shown in Section 11 of [13], this metric g is 
Lorentzian on S x oo). 

Letting p be any point in S", and choosing an open subset G 
which contains p and has coordinate representation {y^, ■ ■ ■ we define 

coordinates {y^, y^, - ■ ■ , y^} on R x S* with = t, and we consider the PDE 
system 

V"Va0 - V'{(l)) = 0, V^F^. = V^F^. + d,{g- %A^) = 0, (209) 

R^.u = d^(t)d,(t) + ■^^V{(t))g^, + F^^F, " - -^-l—^g^.Fp^FP'', (210) 
where (as in Subsection 12. ip 

R^,:=R^, + V(;D,^, V,:=T,-T^, T,:=g,,g^^T''^p, (211) 

G ■= VM^, F^, := F^. + g^,{g - V^A^), (212) 

and r[^^ is the Christoffel symbol corresponding to the reference metric g. 

We may choose gauges so that = Q = Q on S. Since the constraints 
f|T2|) - f|T^ are satisfied on 5, we necessarily have dtV^ = dtQ = on 5*. 
It then follows from local existence and uniqueness results and arguments 
similar to those of [TS', that we obtain a patch of spacetime development 
{Wp, gp, (pp, Ap) which satisfies ([I]) - ([3]); here Wp is a spacetime neighborhood 
containing the point p. 

Step 3: Patching together the spacetime development patches. 
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With a bit of described in Section 16 of [13], it is straightforward 

to combine the development of the data related to U (from Step 1), with the 
developments of data on sets Wp, for a sufficient collection of points p (as 
per Step 2), and thereby construct a Cauchy development (M, g, cj), A) of the 
initial data on E". This completes the proofs of Theorem [1] and Theorem [2j 



8 Conclusion 

In this paper we have proven the asymptotic GFC-stability of the inflationary 
cosmological model 



X 2 Co 

as a solution to the Einstein-Maxwell-Scalar field equations 

1 



R 



T 

-J- II 



V^V^^ - V'{<P) = 0, 



0, 



with 



P 



(n- 1)A2' 



Co = In 



2{np 



9.0 + y(0)) + (F^.F, 



There are other model solutions of similar field equation systems for which 
one might be able to obtain similar conclusions. These include spacetime so- 
lutions of field equations including charged scalar fields or multiple scalar 
fields. Along with considering various fields coupled to the Einstein equa- 
tions, one might also consider the spacetimes with non-fiat spatial slices. 
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